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ABSTRACT 

The evolution of radiation emitted during the dynamical collapse of metal- 
free protostellar clouds is investigated within a spherically symmetric hydro- 



(N 
> 

in 

dynamical scheme that includes the transfer of radiation and the chemistry of 
C^i ' the primordial gas. The cloud centre collapses on a time scale of ~ 10 years, 



thanks to line cooling from molecular hydrogen (H2). For most of the collapse 
<"^ ■ time, when the evolution proceeds self-similarly, the luminosity slowly rises 

up to ~ 10 36 erg s _1 and is essentially due to H2 IR line emission. Later, con- 



tinuum IR radiation provides an additional contribution, which is mostly due 

^ to the accretion of an infalling envelope upon a small hydrostatic protostellar 

core which develops in the centre. We follow the beginning of the accretion 
• <-^ ■ 

^^ ' phase, when the enormous accretion rate (~ 0.1 Afoyr" 1 ) produces a very 

high continuum luminosity of ~ 10 36 erg s _1 . Despite the high luminosities, 
the radiation field is unable to affect the gas dynamics during the collapse 
and the first phases of accretion, because the opacity of the infalling gas is 
too small; this is quite different from present-day star formation. We also find 
that the protostellar evolution is similar among clouds with different initial 
configurations, including those resulting from 3D cosmological simulations of 
primordial objects; in particular, the shape of the molecular spectra is quite 
universal. Finally, we briefly discuss the detectability of this pristine cosmic 
star formation activity. 

Key words: stars: formation - stars: mass function - cosmology: theory - 
galaxy: formation 
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2 E. Ripamonti, F. Haardt, A. Ferrara, & M. Colpi 
1 INTRODUCTION 

As the temperature of the cosmic bath decreases, atoms start to recombine and therefore 
decouple from the CMB radiation at redshift ~ 1100. The baryonic Jeans mass immediately 
after this event is given by 

Mj ~ 4.5 x 10 4 (5*) (Qh 2 )-^M Q , (1) 

where Q and 0& are the matter and baryon density parameters, respectively. Masses larger 
than Mj are gravitationally unstable and should, in principle, collapse. However, in order 
for the actual collapse to occur a more severe condition must be satisfied, i.e. that the 
cooling time of the gas is shorter than the Hubble time at that epoch. In fact, radiative 
losses provide the only way for the gas to lose pressure and settle down in the potential 
well of the host dark matter halo. Since the virial temperature corresponding to the masses 
of the first objects (PopIII) is typically £ 8000 K, cooling by hydrogen Lya excitation is 
strongly quenched, and the only viable coolant in a primordial H-He plasma is molecular 
hydrogen. H 2 is produced during the recombination phase, but its relic abundance is very 
small (/h 2 ~ 10~ 6 ). This primordial fraction is usually lower than the one required to 
fulfill the above time-scale constraint, but during the collapse phase the molecular hydrogen 
abundance reaches values that are high enough to allow the collapse to continue. It follows 
that the fate of a virialized lump depends crucially on its ability to rapidly increase its 
H 2 content during the collapse phase. 

Haiman, Thoul & Loeb (1996) and Tegmark et al. (1997) have addressed this question by 
calculating the evolution of the H 2 abundance for different halo masses and initial conditions 
for a standard CDM cosmology. They conclude that if the prevailing conditions are such that 
a molecular hydrogen fraction of order of /h 2 ~ 5 x 10~ 4 can be produced during the collapse, 
then the lump will cool, fragment and eventually form stars. This criterion is met only by 
the largest PopIIIs, implying that at each virialization redshift one can define a critical mass, 
M cr i t) such that only protogalaxies with total mass M > M cr i t will be able to eventually 
form stars. As an example, a 3a fluctuation of a Cold Dark Matter primordial spectrum, 
has M cr u ~ 10 6 M Q and collapses at z ~ 30. Other studies have revised the value of M cr it 
as a result of a more accurate treatment of the micro-physics (Nishi & Susa 1999; Fuller & 
Couchman 2000; Machacek, Bryan & Abel 2001); however, these results do not qualitatively 
alter the basic argument given above. 
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Radiation from the first forming stars 3 

Recently, some works have started to tackle the formation and collapse of PopIII objects 
through numerical simulations (Abel et al. 1998; Nakamura & Umemura 1999; Nakamura 
& Umemura 2001; Bromm, Coppi & Larson 1999, 2001; Abel, Bryan & Norman 2000 - 
ABN00; Bromm, Coppi & Larson 2001 - BCL01; Bromm et al. 2001) based on hierarchical 
scenarios of structure formation. These studies have shown that gravitational collapse in- 
duces fragmentation of the first pregalactic objects with initial baryonic mass ~ 1O 5 M into 
smaller clumps of typical mass of about 1O 2_3 M , which corresponds to the Jeans mass set 
by molecular hydrogen cooling. 

Tracking the subsequent evolution of these clumps is a very challenging problem as it 
requires the simultaneous solution of the hydrodynamic equations and of radiative transfer 
in lines. At present, due to the tremendous computational demand, this as only been possible 
in one dimensional formulations of the problem. Star formation in the early universe naively 
appears easier to understand as compared to the present, because several complicating effects 
can be neglected to a first approximation: among these are magnetic fields, dust grains and 
metal enrichment. 

In their pioneering work, Omukai & Nishi (1998) (ON98) showed that at a certain stage 
of the collapse a quasi-hydrostatic central core is formed whose specific mass is almost 
independent on the details of the problem, with a typical value of ~ 5 x 1O _3 M . However, 
as the accretion rate of the infalling gas is found to be rather high, around 1O~ 2 M yr _1 , 
the core mass at its formation might be unimportant in setting the final mass of the star. 
In the absence of any effect quenching accretion, a large fraction of the initial object can 
become soon part of the protostar. Pushing this conclusion a bit further, one might predict 
a top heavy IMF for this first generation of sources. 

If the infall can instead be stopped, different scenarios can be envisaged. It has yet to be 
firmly established if standard ways proposed to halt the infall continue to work under pri- 
mordial conditions: radiation force might be opacity limited; bipolar flows need some MHD 
acceleration process and therefore seem to be excluded by the weak magnetic field present 
at the main formation epoch (Gnedin, Ferrara & Zweibel 2000), although the magnetic field 
could be amplified during the protostar formation process. Angular momentum barrier re- 
lated to disk formation and competitive accretion among different clumps are slightly more 
promising, but much more difficult to model accurately. The estimate of the first mechanism 
(radiation reaction) requires to accurately derive the properties of the radiation field pro- 
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duced during the collapse and to evaluate the effects of the radiation force on the infalling 
envelope. This is one of the main aims of the present paper. 

The second motivation for our study consists in assessing the detection chances of distant 
PopIII objects via the (infrared) radiation they release during the collapse phase. Ciardi 
& Ferrara (2001) have already pointed out that IR H 2 molecular lines produced by the 
cooling of supernova shocked PopIII gas will be very likely observed by the Next Generation 
Space Telescope. However, as H 2 roto- vibrational lines are essentially the only carrier of the 
entire gravitational energy release, it seems worthwhile to compute the intensity, shape and 
evolution of the emitted spectrum in detail. Recently, Kamaya & Silk (2001) (KS01) have 
estimated the H 2 luminosity from primordial collapsing stars using a very approximated 
treatment. 

To investigate the described problem we have performed ID numerical simulations in- 
cluding a detailed treatment of the molecular line radiative transfer. The study is germane to 
the one of ON98 but it includes the radiative force term in the dynamical equations adding 
a number of imporvement, as we outline throughout the paper. 

The paper is organized as follows. In §2, we present our numerical methodology. The 
choice and justification for the initial conditions, and the results of the simulation are dis- 
cussed in §3. In §4 we discuss the luminosities and the detect ability, and in the final section 
we give a brief summary. 

2 NUMERICAL METHODS: THE CODE 

We developed a 1-D Lagrangian hydrodynamical code, including full treatment of radiative 
transfer and chemical evolution (with careful treatment of the H2 molecule, which is the 
main coolant in the initial phases of the collapse). Such code is similar to the one described 
by ON98. Here we only summarize the relevant differences between the two; a more detailed 
description is given in the Appendix. 

2.1 Dynamics: radiative force 

We included the radiative force term into the momentum equation, which then becomes 
dv n dp GM,, 

— = — 47IT — — : 

dt dM r r 



, ap UlVlr „ ,^s 

- 4 ^ 2 ^iV " -ZJT + f«*> ( 2 ) 



where we are exploiting the spherical symmetry, and where 
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Radiation from the first forming stars 5 

/rad = - / K v F v dv (3) 

is the radiative force per unit mass, F v is the specific energy flux at frequency v (both / ra d 
and F v are considered in the radial outward direction), and k v is the opacity at the frequency 

v. 



2.2 Chemistry and thermodynamics: non ideal effects 

As discussed by several authors (e.g. Hummer & Mihalas, 1988; Mihalas, Dappen & Hum- 
mer, 1988; Saumon, Chabrier & Van Horn, 1995; hereafter SCVH95), at high densities 
(p £ 10~ 3 gem -3 ) particles are jammed closely together, and the bound electron orbitals 
filling too large a volume fail to survive so that electrons migrate from atom to atom. Such 
an effect is known as pressure ionization and generally it is not taken into account in the 
Saha equations. As an example, a pure hydrogen gas with p = 1 g cm -3 and T = 3 x 10 4 
K is largely ionized (the mean distance between nuclei is only twice the Bohr radius of an 
H atom on the fundamental level, and the electrons are free in a Fermi gas, interacting 
simultaneously with several ions), while solving the Saha equations results in a completely 
molecular gas, which is clearly unphysical. In such regimes, the gas is also non-ideal, and 
the perfect gas equations are poor approximations. 

For these reasons, at high densities and temperatures we find the chemical abundances 
and thermodynamical properties (such as the pressure p, the internal energy per unit mass 
u and the entropy per unit mass s) by interpolation of the values tabulated by SCVH95, 
which are calculated using the free energy minimization technique. SCVH95 calculations 
include several non-ideal effects, such as, (a) the partial electron degeneracy, (b) the Coulomb 
interactions among charged particles, (c) the finite extension of "hard-sphere" particles, 
(d) weak attractive interactions (e.g. van der Waals), and (e) modifications of the internal 
partition functions due to the shift of the energy levels. 

3 THE SIMULATIONS 

3.1 Initial conditions 

The main parameters specifying the initial condition of our runs are: the total cloud mass 
M to t, the temperature, density, and velocity profiles T (M r ), p (M r ), t> (M r ), respectively, 
and the hydrogen molecular fraction /h 2 . 
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Figure 1. Spatial structure of the cloud in run P100/2 at 10 different stages of collapse (see Table El). Panels refer to density 
(a), infall velocity (b), temperature (c) and to the fraction of H nuclei inside H2 molecules (d). On the abscissa we use the 
mass enclosed within radius r (M r ). Shell edges are denoted by small tags onto the curve referring to the initial configuration; 
in panel (d) we omit stages 6 and 8 for clarity. The small oscillations in the temperature profile of stage 9 around M r ~ — 2 
are due to numerical effects. 

Several previous investigations set constraints on such parameters. Theoretical predic- 
tions of the mass of fragments leading to the birth of the first stars are quite uncertain, 
ranging from a minimum of ~ 0.1 M (Palla, Salpeter & Stahler 1983), to £ 200 M 
(Hutchins 1976). On the other hand, there is agreement about the value of /h 2 ~ 5 x 10~ 4 , 
and about the value of T (M r ) which is of order of few hundreds degrees (e.g. Tegmark et 
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Table 1. Initial conditions 



Run 


P0,c 


To, c 


/h 2 


Notes 




(mjf cm -3 ) 


(K) 


(10- 4 ) 




P100 


2.6 X 10 6 


270 


5 


cfr. ON98(A,B) 


P100/2 


2.6 x 10 6 


270 


5 


refinement of P100 (140 shells) 


P1000 


2.3 x 10 4 


246 


5 




P10 


3.5 x 10 8 


249 


50 


cfr. ON98(C) 


1100 


7.8 x 10 7 


300 


5 




11000 


1.1 x 10 7 


300 


5 




N1500 


2.0 x 10 6 


440 


5 


from ABNOO 



The letter in the label of each run gives the assumed initial temperature/density profile: polytropic, isothermal (I) , and numerical 
(N), i.e. after ABNOO results. It is followed by a number which gives the assumed total cloud mass Mtot, in solar masses. 

al. 1997). Detailed numerical simulations (BCL01, ABNOO) suggest that the runaway col- 
lapse starts soon after the (number) density reaches a value of 10 4 — 10 5 cm -3 . Such value 
is determined by the transition from non-LTE cooling (with A oc n 2 ) to LTE cooling (with 
A « n), slowing down the contraction, and allowing the gas to reach a quasi-hydrostatic 
state. This transition occurs when T ~ 300 — 400 K and /h 2 ~ 5 — 10 x 10 -4 , inside a cloud 
with a mass of a few thousand M & , a value close to the Jeans mass of the cloud. 

To check how our results depend on the choice of the initial conditions, we run a set of 
models with different initial mass M tot , ranging from 10 to ~ 1500 M & . We use a constant 
profile for the H2 abundance and zero velocity (the cloud is at rest). The density and temper- 
ature profiles inside the cloud are derived imposing hydrostatic equilibrium and a polytropic 
index equal to 5/3 or 1 (isothermal). In addition we have considered 'numerical' profiles as 
given by ABNOO from cosmological simulations. In runs P10, due to the high value of the 
initial central density, we used a higher value for /h 2 • The initial electron fraction adopted is 
f e - = 10~ 10 . With the exception of run P100/2, we always divide our simulated cloud into 
100 Lagrangian shells, with massive external shells and much lighter internal ones (e.g. , in 
run P100 each external shell comprises ^ 1 M© of gas, while the central one has a mass of 
only 3 x 10 -7 M©). Such choice is aimed to follow with good resolution the fast evolution of 
the centre of the cloud, which accounts only for a tiny fraction of the total mass. Details of 
the initial conditions and general set up are given in Table [|. Note that runs P100 and P10 
are directly comparable to those of ON98. Run P100/2 is identical to run P100, but the ob- 
ject is divided into 140 (rather than 100) shells; the "extra" shells were used to improve the 
resolution in the particularly important region where the shock forms (—3 ^ logM r £ — 1). 
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Figure 2. Evolution of density and temperature of the innermost shell (thick line). Points A and B mark conditions when 
the cloud becomes optically thick and when the hydrostatic core forms, respectively. Thin dotted lines show the loci of points 
where equilibrium H2 fraction is 0.9 (bottom), 0.1 (middle) and 0.01 (top). Dashed line show the theoretical locus of points 
where a generical cloud should become optically thick. 

3.2 Dynamical evolution 

The dynamics of the collapse is illustrated in Figure [I] that shows the evolution of the cloud 
in run P100/2 at 10 different stages (see Table §). In Figures and § we show the evolution 
of the temperature and of the main hydrogen species in the most central shell as a function 
of density. In the rest of this section we will describe the results obtained in run P100/2; 
results of other runs are quite similar. In runs P100, P100/2 and P10 we find good agreement 
with the results of ON98. 

In the beginning, the evolution is very well described by a Larson-Penston self-similar 
solution (Larson 1969 - L69; Penston 1969) with a density profile scaling as p oc r~ a (with 
a ~ 2.2). The self-similar evolution breaks down between stages 6 and 7, when a quasi- 
hydrostatic core forms in the centre and starts growing in mass. 



3.2.1 The collapse phase 

During collapse the central regions of the cloud undergo several important chemical and 
thermal transformations. When the number density reaches £ 10 8 cm -3 (stage 1), 3-body 
reactions become efficient, converting the bulk of H into H2 (Fig. p], stages 2-3; Fig. [3|). 
The growing H2 fraction increases the gas cooling rate, only partially counterbalanced by 
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Table 2. Collapse stages (run P100) 



Stage 


T c 


At 


(</ - U) 


Notes 




(K) 


(yr) 


(yr) 







270 


— 


740000 


Initial conditions 


1 


150 


720000 


20000 


Start Rad. Trans. 


2 


650 


19000 


930 




3 


1500 


920 


9.6 




1 


3000 


9.2 


0.40 


H2 is dissociating 


5 


6500 


0.060 


0.34 




6 


15000 


0.005 


0.34 


Core is forming 


7 


30000 


0.002 


0.34 


Core has formed 


8 


65000 


0.011 


0.33 




9 


100000 


0.33 


— 


Final stage 



radiative transfer effects (see next section 3.3). The evolution of the thermal properties of 
the gas is shown in Figure f|, where the adiabatic heating T a d, the radiative cooling rate A rac j 
and the effective heating r eff = T ad + r ch (r ch is the heating - or cooling - due to chemical 
reactions) are plotted vs. the gas temperature. All these quantities refer to the central 
shell. As T ;$ 1550 K, the chemical heating due to the H2 formation dominates over r at j 
and heating and cooling terms almost balance, leading to a very slow temperature growth. 
For T ^ 1550 K H2 is destroyed rather than formed, and the chemical term becomes a net 
cooling (this is the reason why the curves r e g- and r a d cross at T ~ 1550 K). At this stage the 
chemical and radiative cooling terms are of the same order. Radiative cooling is dominated 
by line radiation for T ;$ 1800 K, while continuum emission is important above. However 
continuum radiative cooling soon becomes highly inefficient because of the large optical 
depth: for T £ 2350 K the innermost regions become almost adiabatic, and compressional 
heating is only matched by chemical cooling. 

We just mentioned that H2 starts to dissociate when T ^ 1550 K. This process is slow 
because of the density increase and the chemical cooling associated to the dissociation process 
(4.48 eV per dissociated molecule are lost by the gas thermal energy). As a result, the 
H 2 fraction at the end of the self-similar phase (T c ~ 15000 K, stage 6) is still ~ 5%, even 
in the hottest central regions. 

3.2.2 Core formation 

During the last stages of the self-similar collapse (stages 5-6; evolution from point A to point 
B in Fig. |2|), when the centre of the cloud is essentially adiabatic, H 2 dissociation acts as a 
"thermostate" , preventing a faster increase in temperature and pressure, which would halt 
the collapse. In fact, dissociation absorbs ^ 50% of the heat generated by gas compression, 
which (due to the opaqueness of the cloud's centre) could not have been dissipated otherwise. 
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Figure 3. Evolution of the fractions of the main hydrogen species in the innermost shell as a function of density: H2 (solid 
line), H° (dashed line) and H + (dotted line). 
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Figure 4. Radiative cooling rate (dashed), adiabatic heating rate (dotted) and effective (i.e. adiabatic plus chemical) heating 
rate (solid) in the innermost shell as a function of its temperature. The left vertical line marks the transition from H2 formation 
to H2 dissociation, while the right vertical line separates the optically thin and optically thick regimes. 

When nearly all H 2 molecules are dissociated, the thermostatic effect is no longer effective, 
and the central temperature and pressure rise rapidly stopping the collapse of the central 
regions. At the end of H 2 dissociation (T c ~ 20000 K) the chemical energy and the thermal 



© 0000 RAS, MNRAS 000, 000-000 



Radiation from the first forming stars 11 



£.-2 



BO 
O 



- 

H + 

; 


\ / 


,i-'-j 


— 


1 


A 

/ \ 

/ \ 
/ \ 
/ \ 




- 


H° v "' 


- 


i 

i 




: 


■ 1 


; \ 




^H 2 


- 


- \ 


, , i 


,,l, 


- 



-6 -4-2 2 

Log(M r ) [MJ 

Figure 5. Chemical structure of the cloud at stage 9: H2 (solid line), H (dashed line) and H + (dot-dashed line); the ionized 
hydrostatic core, the inner molecular and the outer atomic envelope are quite apparent; the peak of H° curve at logM r ~ — 1.3 
is located in the region where the (smoothed) shock brings down the infall velocity of the shell from the maximum value of 
~ 15 kms -1 to ~ 2 kms - 1 . The small oscillations for —2.5 ^S logM r ^S — 1.5 are numerical. 

energy per H° atom are similar (~ 2.2 eV and ~ 2.6 eV, respectively). This fact reflects the 
virial theorem: half of the cloud gravitational energy is converted into thermal energy after 
compression, while the rest is tapped into chemical energy through dissociation. 

Outside the core, the gas is still collapsing with a velocity close to free-fall. This results in 
a shock front forming at the edge of the core, where the infalling material is suddenly halted 
(see Fig. |I].b). Unfortunately, the numerical method (i.e. the introduction of an artificial 
viscosity) used for the treatment of hydrodynamical effects is inaccurate when dealing with 
strong shocks, since they are spread upon several shells. However, though the treatment of 
the shock physics is only approximate, the general picture of the cloud evolution is essentially 
correct. In run P100/2 we also improve the situation by refining the shells in the region where 
the shock forms. 

The core evolution is driven by two main processes. First, as can be seen in Figures 
|l|.b, and || the core mass (defined as the mass inside the radius where the infall velocity is 



< O.lfff, vg = J2GM r /r) grows, due to the accretion of infalling material. Second, even if 
the infall velocity is drastically reduced, the core is compressed by the infalling material, 
and then heats up (as can be seen in Fig. |2|, beyond point B). 
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3.2.3 Cloud structure after core formation 

After core formation and during its subsequent evolution, the H2 fraction in the centre never 
rises again to important levels^. H2 tends to form a layer of several solar masses between 
the small core (where H is progressively ionized) and the outer layers (where H is essentially 
atomic, see Figure |5|). After stage 6 we can recognize three different regions, i.e., the core, 
the molecular envelope, and the outer region (see Figure ||]). 

(1) The central core is a region whose infall velocity is much smaller than the free-fall 
value. At core formation, the gas here is basically neutral, but at stage 9 it has become 
mostly ionized (but nu+/nn is still only ~ 2) because of the pressure ionization process. 
The outer edge of the core is marked by a shock front, where material accretes nearly with 
free-fall velocity. Since the artificial viscosity technique "smooths" the shock, we have a 
thin transition region, where neutral H is the dominant species; that may be an numerical 
artefact. 

(2) Outside the shock defining the core outer edge, there is an envelope of ~ 3 M composed 
by molecular hydrogen (the shock also marks a sharp transition in chemical composition, 
from H2 to H° or H + ). Such envelope is essentially in free fall (v/vg is in the range 0.3-0.7, 
decreasing outward). Its outer edge is roughly at a radius ^ 10 15 cm. 

(3) The remaining outer gas, accounting for most of the cloud mass, has a chemical com- 
position only slightly different from the initial one, since reactions leading to H 2 formation 
are slow (mainly because of the low density), and involve only a minor fraction of H atoms. 
This layer is the extension to the inner molecular envelope, and has a velocity v /vg between 
0.1 and 0.3. 

At stage 9 the core temperature and density are ~ 10 5 K and ~ 0.3 g cm -3 , respectively. 
In these conditions the Courant time-scale is extremely short ( ^ 300 sec), and it becomes 
more and more difficult (and costly, in terms of CPU time) to follow the evolution further. 
For this reason, we stopped our calculations when the central temperature reached ~ 10 5 K. 



The central value of /jj 2 at late stages is quite uncertain: using SCVH95 results we obtain /u 2 ~ 0.03; had we used results 
such as those of Mihalas, Diippen & Hummer 1988, we would have obtained /jj < 0.01; ON98 find /jj ~ 0.4, but this is likely 
to be too high since they do not consider pressure ionization effects. 
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3.2.4 The accretion phase 

After core formation, the evolution continues because of the accretion of the outer envelope 
onto the core. 

Despite our approximate treatment of the shock, we find that the accretion rate, as in- 
ferred from the full calculation, is extremely high: the core grows from ~ 0.3 x 10~ 2 M & when 
it forms, to ~ 4.5 x 10~ 2 M & in ~ 0.30 yr, which means M ~ 0.14 M yr _1 . Such accretion 
rate is much larger than typical estimates for present-day star formation (M £ 10~ 4 M Q 

In order to extend our results to a more advanced state of accretion, we modified our 
code "freezing" the evolution of the shells deeply embedded in the core. In practice, when 
a shell has i) a small infall velocity (v/v$ < 10 -3 ), ii) a high temperature (> 5 x 10 4 K), 
and iii) is located at least 5 shells inside the shell where the infall velocity is maximum, 
we completely stop its evolution (i.e., we set velocity at and keep all the other quantities 
constant to the last computed values). 

This greatly speeds-up our calculations because there is no more need to follow the shells 
with the shortest time-scales, drastically reducing the number of integration steps needed. 
Such "freezing" is appropriate, since these shells do evolve on time-scales much longer than 
those imposed by the Courant condition and have little influence upon the layers where 
accretion occurs. There are two major concerns: 

- the neglect of the secular changes in the core radius: for example, in the full calculation 
the radius of the 10 innermost shells undergoes a reduction of about 40% (from 1.6 x 10 10 
cm to 1.0 x 10 10 cm) as T c goes from 50000 K (stage 8) to ~ 10 5 K (stage 9), while in the 
simplified calculation performed between stages 8 and 9, we keep this radius constant. This 
change is less important in more external shells: the difference in the shock position between 
the "full" and the "frozen" run at stage 9 is found to be ~4%. 

- the radiation produced upon the shock surface could ionize a thin shell outside the shock 
itself, as reported in Staher, Palla & Salpeter (1986) (SPS86), modifying the dynamics of 
accretion (see next sections). 

Instead, the influence of the lowered temperatures in central shells is completely negligi- 
ble, since the optical depth from the "frozen" shells is very high (r c £ 10 12 ). 

Figure || shows the result of the accretion calculations performed with the modified 
code: we stop ^ 10 years after core formation. In the first 0.34 years this results are almost 
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Figure 6. Evolution of the core mass. The three lines with dots refer to model P1000 (top line), P100 (middle line) and P10 
(bottom line). The dashed line shows the predictions of ON98 (eq. 11), while the dotted one shows the results of the calculation 
described in the text, using the parameters of run P100. 

identical to those obtained with the "complete" code; instead, points at later times should 
be considered indicative. 

ON98 investigated the evolution of the mass accretion rate on the newly formed hy- 
drostatic core using the Larson-Penston self-similar solution (L69, Penston 1969) and its 
extensions after core formation (Hunter 1977, Yahil 1983, Suto & Silk 1988). For a model 
close to our model P100 (and P100/2), they found that the core mass grows following a 
power-law 

M.(t) = 0.11 M ((t - t cf )/yr) a73 . (4) 

We also find an empirical model which could better fit the numerical results, based on 
the fact that in our calculations we find that the mass flux Mf{ nx (M r ) = 47rr 2 /w m is well 



approximated by a power-law. Such empirical model leads us to predict that 

\ 0.657 



M*(t) = 9.2 x Kr 2 Mo I 0.0003 



,yr 



(5) 



Both these analytical expressions are compared with the numerical results in Figure [| 



3.2.5 Importance of initial conditions 

The above description essentially applies to all the models we computed: the evolutionary 
sequence is always the same, even if there are small quantitative differences. Models lose 
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Figure 7. Comparison of the evolution of central conditions (density and temperature) in the several of our computed models: 
from top to bottom, the lines refer to model P1000, 11000, P100, 1100 and P10; model N1500 is almost indistinguishable from 
model P100, except in the first stages. As in Fig. H, points A and B denote the conditions where (in model P100) r first becomes 
> 1 and where core formation occurs, respectively. 

memory of the initial conditions almost completely during the self-similar phase, as can be 
seen in Fig. [F], where tracks referring to the different models are very close for n c > 10 14 cm~ 3 , 
i. e. well before the cloud becomes optically thick (point A) or forms the central hydrostatic 
core (point B). 

There remain only relatively small differences (within a factor of 2 or 3), in particular: 

- Models with higher total mass (M tot ), have higher temperatures than low M to t models, 
for any given central density n c ; 

- Models with high M tot have cores that grow faster than low M tot models (see Fig. ^). 

Both these trends can be easily explained: when we impose that in the initial stage 
the clouds are close to equilibrium, we implicitly force high-mass models to have higher 
temperatures at any given density; even if quantitative differences are strongly reduced in 
the self-similar phase, high-mass models remain warmer than low-mass models at all stages. 
Furthermore, the accretion rate depends on temperature, since M ~ M Jeans (T, p)/t ff (p) oc 
T 3 / 2 (see Stahler, Shu & Taam, 1980) and temperature differences explain also those in M. 

The behaviour of the N1500 model (the only one for which we do not impose quasi- 
equilibrium initial conditions) is particularly relevant: after some "wandering" in the n c — T c 
plane, its central conditions evolve almost exactly along the evolutionary track of model 
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Figure 8. Radial outward optical depth as a function of enclosed mass for (a) continuum radiation and (b) for radiation at 
the center of the roto-vibrational H2 (0,2)— >(0,0) transition (A = 28.3 /im); the top curve in this panel, taken when r c = 1, is 
represents all stages between 4 and 9. 

P100, indicating that the collapse does not preserve the memory of initial conditions and 
the cloud evolves along an equilibrium track. 

3.3 The radiation field 

The evolution of the radiation field is intimately coupled with dynamics and our physical 
description of forming primordial stars fully accounts for their emission properties. Here we 
give details on the evolution of the optical depth, luminosity and spectra during the two 
phases of collapse and core formation, and at the very beginning of the accretion phase. 



3. 3. 1 Luminosity from the collapsing cloud 

At the onset of dynamical evolution, radiation results from the cooling of the optically 
thin H 2 present in the cloud. At stage 1, as 3-body reactions keep converting H atoms into 
H 2 molecules, the line optical depth rapidly grows, reaching significant values. At stage 2 
there are already 14 lines with r > 0.1 (although none of them with r — 1), and this 
number rises to 42 lines (half with r > 1) by the time the central temperature reaches 1000 
K (between stage 2 and stage 3). In Figure ||a we show, as an illustrative example, the 
evolution of the optical depth profile for a selected line during the self-similar phase. 

The increase of the optical depth obviously affects the gas cooling properties (as can be 
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seen from the flattening of the radiative cooling curve in Fig. ^), but it does not influence the 
dynamics in a fundamental way. The increase in H2 abundance is associated with an increase 
of the gas temperature and of the number of relevant transitions of the H2 molecule. As the 
temperature grows, new, high quantum number lines are initially optically thin, providing 
a minimum level of cooling. Though line cooling is severely quenched by opacity (e.g., at 
T c ~ 2000 K only a fraction ~ 0.02 of the line flux is not re-absorbed locally) its absolute 
value increases anyhow as T c increases; this prevents the temperature from rising too fast, 
possibly stopping the collapse in its early phases. 

As already mentioned, when T c = 1000 K we start to include cooling from the continuum, 
which rapidly gains importance, becoming the dominant process in the central ~ 0.1 M Q 
already when T c ~ 1800 K. The evolution of the continuum opacity profile is shown in Figure 
|8|b, where the formation of the core sharp edge is apparent. 

The corresponding emitted luminosities vs. time are plotted in Figure ^. Note that even 
when continuum radiation dominates the central emission, the H2 lines still have a compa- 
rable or slightly larger overall luminosity. This is due to the fact that line radiation comes 
from the outer, massive H 2 envelope, while continuum comes from an inner, much smaller 
region (only slightly larger than the core itself); this can be seen in fig. (TO, where we show 



the radial dependence of the emitted luminosity; in particular, note that the drop in the 
central line luminosity has negligible influence upon the total line luminosity, and that the 
growth of continuum luminosity stops at M r ~ 10~ 2 — lO^ 1 M , while line luminosity is still 
growing at M r ~ 1 M Q . 

3.3.2 Spectra 

In Fig. [O], we show the emitted spectra at 3 different stages. In the upper panel the spectrum 
comes from line emission only, while in the lower panels we assumed that the continuum 
spectrum is a black-body at the effective temperature of the central shell (1500 K - panel 
b) and of the continuum photosphere {i.e. of the region where the radial outward optical 
depth is ~ 1). Such temperature remains roughly constant at ~ 2300 — 2400 K from the 
birth of the continuum photosphere to the end of our calculations. 

In appendix B we give a full listing of the luminosities of all the calculated lines, plus 
the integrated continuum luminosity at 4 stages (the 3 depicted in fig. |TT| plus the stage 
when T c = 1000 K). We note that the durations of the various stages are very different: line 
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luminosities resemble those listed in the Lq^q column (corresponding to fig. p[a) for ~ 10 4 yr, 



while they resemble those given in columns L WQ0 , L1500 and L Tc= i only for ~ 300 yr, ~ 40 yr 
and ~ 2 yr, respectively!]]. 

We note that the spectra are produced only in the central collapsed region: line emission 
comes essentially from the molecular envelope, while continuum emission is produced close 
to the surface of the hydrostatic core. In practice, the whole protostellar luminosity comes 
from the central ~ 3M . The structure of this central region is quite similar in all the models, 
since it is largely determined by the self-similar collapse, rather than by initial conditions. 
For this reason, the spectra emitted by different clouds during these phases should be similar. 

3. 3. 3 Radiation pressure effects 

We have found that the role of radiative force is negligible during all the evolution we 
followed, e.g. at least until the first phases of accretion (our stage 9). In a test run we 
have removed the radiative term in the momentum equation and re-calculated evolution. 
We then compared the results with those of the "full" code, and found that the difference 
is completely negligible. 

Returning to the early phase of accretion, despite the very high luminosity (~ 10 36 erg s _1 
in the continuum, plus other ~ 10 36 erg s _1 in H 2 lines), which is comparable to the Ed- 
dington luminosity L Edd — 1-25 x 1O 36 [M/(O.O1M )] erg s _1 , the opacity of the infalling gas 
is much smaller than the electron scattering opacity (k cs = er T /mH — 0.4 cm 2 g _1 ), since 
the infalling material is essentially molecular, with a very low ionized fraction, and other 
opacity sources (Rayleigh scattering, bound-free and free-free transitions, collision-induced 
absorption of H2 ) give small contributions. We find that the opacity is everywhere less than 
0.05k os at all times before the formation of the hydrostatic core; typical values are actu- 
ally several order of magnitude less than this upper limits. The same is true also after the 
formation of the hydrostatic core, if we consider only the accreting layers. 



t The last duration (~ 2 yr) is actually only a lower limit, since line emission is still at high levels when our calculations stop; 
however, the H2 envelope where line emission originates is not likely to last for ^ 100 yr after core formation, since it could 
be accreted onto the core, or photodissociated because of the outgoing radiation. 
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3.3.4 Bridging collapse models with accretion models 

SPS86 and recently Omukai & Palla (2001) studied the subsequent phase of accretion by 
modeling the hydrostatic core as a star of increasing mass, surrounded by an infalling en- 
velope whose structure was determined by an assumed steady mass accretion rate M* = 
4.4 x 10~ 3 M yr _1 . They found that the shock surface is inside an opticaly thick ionized 
region, whose outer edge is roughly coincident with the photosphere. They also found that 
between the shock front and the photosphere the radiative force term is important, since 
K ~ k cs and the luminosity is not far from the the Eddington value. Despite this fact, both 
works find that the radiative force is not able to stop the accretion before the core mass 
grows to very high values (> 100 M ). 

We find that also in our models the shock front is inside the optically thick region, but 
the gas between the shock and the photosphere is not ionized, mainly because the density 
in the region is relatively high (~ 10~ 6 g cm -3 ), about 1000 times larger than the value 
that can be estimated from Fig. 12 of SPS86; instead, we find that the temperature is only 
slightly lower than in the SPS86 model. The ionization fraction inside such region is low and 
remains such ( ^ 0.05) even if we make the hypotesis that radiation coming from the shock 
heats the gas inside the photospheric shells, so that gas and radiation temperature are the 
same (~ 12000 K, the post-shock temperature found both by ours and SPS86 calculations). 
So, this gas still has a low opacity, and the radiative force is small. 

Although numerical effects could play some role in this discrepancy, the most likely 
explanation is that the models of SPS86 are admittedly indicative for core masses ^0.1 M , 
i. e. during the "decay of transients" phase, when the arbitrarily chosen initial conditions 
are still important; instead in our models the core mass is always in this range (here we do 
not consider the simplified calculations aimed to the study the growth of the core mass). It is 
quite possible that both results are correct and that a transition occurs after the last stage of 
our calculations. For instance, we note that the difference in the density of the photospheric 
layers (which keeps the ionization fraction at low levels) is mainly due to the large difference 
(a factor ~ 100) between the accretion rate we find and that assumed by SPS86; since in our 
models we also find that the mass accretion rate is decreasing, such difference will reduce. 



© 0000 RAS, MNRAS 000, 000-000 



20 E. Ripamonti, F. Haardt, A. Ferrara, & M. Colpi 



00 

o 







1 i i i 1 i i i 1 


ii i | i 
total 






X / 


,i^— — 


36 


IK. 

* -^, lines 






/ / continuum 






/ / 


- 






/ / 
/ — i— 


(1.3)-»(0.1) - 


34 




/ / 


- 






/ / / 

/ / 

/ 
/ 


- 






/V i 

1 / 1 


(0.2)-(0.0) " 


32 




- 






!./ i 
..--7 / 
/ ?' 
/i , , , fi , , , i 


, , 1 , 



4 2 0-2 

Log(to-t) [yr] 

Figure 9. Evolution of the luminosity in lines (heavy dotted), continuum (heavy dot-dashed) and their sum (thick solid line); 
thin lines refer to two H2 transition: (0,2)— ►(0,0), A = 28.3 /im (dotted) and (1,3)— ►(0,1), A = 2.12 /an (dot-dashed); they do 
not appear when we "freeze" all line luminosities. Note that the x-axis shows the time left before the end of our calculation. 
The oscillations in the continuum luminosity between 0.5 XL log(to — t) ~ — 1 arc connected with core formation and could be 
partially numerical. 







1 | 1 1 1 | 1 1 1 | 1 1 1 | 


1 1 1 1 1 


36 






























/ / 








* / 








/ 
/ 








-- "// 








' / 




34 




/ / // 
• / ft 


— 






/ / // 
/ / // 
• / // 


- 






/ / / / 

' / / / 








/ / / j 
/ / /I 


- 


32 




/ / /I 

/ / / / 

/ / / I 


- 






/ / / / 








' / / 










/ / / 










/ / / 








" 


/ / 




- 


30 


- 


7 / 




- 


28 




1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 


1I1 



-6 -4-2 2 

Log(M r ) [MJ 

Figure 10. Outward total luminosity in H2 lines (solid) and in the continuum (dashed) as a function of enclosed mass at two 
different times: when T c = 2000 K, t c < 0.1 (thin lines) and when T c ~ 2300 K, t c ~ 10 (thick lines). 



© 0000 RAS, MNRAS 000, 000-000 



Radiation from the first forming stars 21 




0.5 l 

Log(A) H 

Figure 11. Emitted spectra when T c is 650 K (a), 1500 K (b) and 2338 K (i.e. just before we "freeze" lines) (c). Line widths 
are actually smaller (A A/ A ,$ 10 — 4 ) than line thickness in this plot, so that overlaps are only apparent. 

3.4 Comparison with metal-rich protostars 

In this section we discuss the similarities and differences between the PopIII (high redshift, 
zero metal) star formation process we addressed and the PopI-PopII (low redshift, solar 
metallicity) star formation. 



3-4-1 Hydrostatic core 

In our work (and also in that of ON98) the birth of the quasi-hydrostatic core is quite 
different to that described by L69 (and more recently by Masunaga, Miyama & Inutsuka, 
1998 and by Masunaga & Inutsuka, 2000) for a present-day protostar because of the very 
different conditions. In fact, those authors found that a molecular core forms first, and only 
later on, H 2 dissociation leads to the formation of a second more internal core, which is now 
ionized and at much higher temperature and density, as can be seen in table |3|. The two 
cores coexist for a short time, while the second core "absorbs" the first one. 

Instead in our models we find a single core, completely different from the first (molecular) 
core of metal-rich models, but closely resembling the second one. 

The reason of both the differences and the similarities between the two cases is essentially 
explained by an argument due to Omukai (2000). He found that the locus of points in the 
n — T plane where r ~ 1, M ~ Mj eans (T, n) and A rac j ~ T a d (cfr. Silk, 1977) is approximately 
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Table 3. Comparison of quasi-hydrostatic cores initial properties 

P100 L69(I) L69(II) 

M/(10 30 g) 6 10 3 

R/(10 n cm) 1 600 9 

T c /K 2 x 10 4 200 2 x 10 4 

p c /(g cm" 3 ) 1 x 10~ 2 2 x 10~ 10 2 x 10~ 2 

f H2 £ 0.03 ~ 1 ~ 

P100: results of our run P100; L69(I): first core of Larson, 1969; L69(II): second core of Larson, 1969. 

described by a line (L) with T oc n 2 / 5 , which is shown in Fig. |2|; on the right of L, a core 
with mass > Mj eans (T, n) is optically thick. 

Our metal free models become optically thick at point A (remarkably close to L), where 
T c is high and H2 dissociation is immediately effective, preventing the formation of a hy- 
drostatic core. Instead, metallic protostars (Z £ 10 Z ) reach L at temperatures too low 
for H 2 dissociation, and the first (molecular) core forms. Such core evolves adiabatically, i.e. 
along a T oc n 7ad_1 line, which remains close to L, because 7 a d ~ 7ad,H 2 — 7/5. In this way 
cores of all metallicities converge towards the vicinity of point A. After point A, the evolu- 
tion of the central regions is essentially determined by the balance between compressional 
heating and H2 dissociation cooling and no more depends on metal content; thus, metal rich 
and metal free models form a core when H2 is exhausted, close to point B. 

The initial core mass ~3x 1CT 3 M Q (only slightly dependent upon initial conditions) 
agrees with Mj eans (T B ,7iB) (where subscript B denotes the values at point B). Such mass 
is also close to the minimum value reached by Mj eans (T c ,n c ) during the cloud evolution 
(see Fig. |T^). This minimum Jeans mass is generally interpreted as the minimum fragment 
mass {i.e. the mass below which further fragmentation is impossible) and we find it depends 
very weakly upon initial conditions and on metallicity. An important corollary is that the 
minimum fragment mass can not be simply found imposing r £ 1, since we have seen that 
such condition, although necessary, is not sufficient for the formation of a stable hydrostatic 
core. 

It is also interesting to note that in our models there is no sign of instability inside the 
hydrostatic core due, for example, to H ionization, even if the core is mostly atomic at stage 
7 and mostlyly ionized at stage 9. 

3.4-2 Accretion rate 

Even if important core properties, such as the initial mass and the evolutionary path of 
central regions in the n—T plane are almost independent of metal content, the mass accretion 
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Figure 12. Evolution of the central Jeans mass (i.e. the Jeans mass calculated assuming p = p c , T = T c ) for several computed 
models (from top to bottom: P1000, 11000, P100, 1100 and P10). Models P100/2 and N1500 are not shown but they are almost 
indistinguishable from model P100. 

rates are completely different, with metal-free objects accreting at a rate £ 10 3 times larger 
than their metal enriched counterparts. For example, in our P100 model core mass reaches 
~ O.3M in ~ 6 years (see fig. |6|) while in the model of Masunaga & Inutsuka (2000) about 
20000 years are required for the accretion of the same mass (see their Fig. 7). 

As already discussed by several authors {e.g. Stahler et al. 1980; SPS86), the reason for 
this difference is found in the initial conditions, in particular in the higher temperature, 
since for an object with M ~ Mj eans we have 

M ~ (k/m H f/ 2 T*l 2 /G ~ 1.8 x 10- 7 T^ 2 M Q yr" 1 . (6) 

If we take a typical initial temperature of our models ~ 300 K we obtain M ~ 10~ 3 M o yr _1 , 
but if we take T ~ 5 K, as appropriate for present-day protostellar clouds, we get a value 
~ 500 times smaller for M. 

We note that @ gives the mean accretion rate for the complete accretion of the pro- 
tostellar cloud onto the core; should we consider accretion of only a fraction of the cloud, 
equation (^) would give higher values of M because the appropriate reference temperature 
would be higher. 
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3.4-3 Radiation feedback 

In the previous sections we showed that radiation pressure is not effective in stopping the 
accretion flow upon a zero-metallicity protostellar core. The case of metal rich protostars is 
different: luminosities are smaller because of smaller accretion rates, but the presence of dust 
brings the opacity at much higher levels. For example, Stahler, Shu & Taam (1981) estimate 
that at the start of accretion radiative deceleration retards the inflow of 42% from the free 
fall value, even if it can not completely stop accretion (see also the models of Masunaga et 
al. 1998 and Masunaga & Inutsuka 2000, which include radiation pressure effects in their 
momentum equation). In a different context, radiative forces upon dust were also found to 
act decisively in halting the accretion inflow upon massive stars (Wolfire & Cassinelli, 1987; 
Larson & Starrfield 1971). 

For sufficiently high density and temperature, ignition of nuclear reactions occurs in the 
core interior. To check if the effect can, in principle, stop the accretion phase, we compare 
the accretion time-scale with the diffusion time-scale of the energy produced. 

As can be seen in Fig. [8], the interior of the core has optical depth r* £ 10 12 . If energy is 
transported radiatively, the time-scale for its diffusion through the core is tdiff — t*R*/c ~ 
10 5 yr[r*/10 12 ][i? ;| ,/(10 11 cm)], where the core radius R* has a value of the order of 5 x 10 11 
cm. We note that, as T c rises, the opacity k inside the core decreases; anyway, k ^ k cs , which 
leads to a strong lower limit r* ^ 10 10 . 

An identical calculation for present-day protostellar cores leads to a time-scale of the 
same order of magnitude. Instead, the accretion times differ enormously. For example, if 
we extrapolate equation (^), we find that our core would take only t 1Q ~ 40 years after 
formation to reach a mass of M* ~ 1M Q (but we note that ti Q ~ 250 yr also in the models 
of SPS86), while in the model by Masunaga & Inutsuka (2000) more than t 10 ~ 1.4 x 10 5 
years are required for the core to reach the same mass. 

So, £i /£diff ~ 1 in metal-rich cores, but ti /tdiff "C 1 in the metal-free case. Accordingly, 
in the metal-free case the energy generated in the core centre (for instance, because of 
ignition of nuclear reactions) cannot affect the accretion, while in the metal-rich case such 
energy could be able to reach the surface and change the conditions of the infalling envelope. 

The arguments discussed in the present section are subject to important uncertainties. 
For example, we did not investigate any energy transport mechanisms, such as convection, 
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different from radiative diffusion, which could be much more efficient (e.g. , see the "lumi- 
nosity wave" described in SPS86). 

4 LINE FLUX AND DETECTABILITY 

The work by KS01 has addresses the important problem of detectability of the emission 
from collapsing primordial protostar. They employed approximate analytical calculations 
(partially derived from the ON98 work) in order to estimate order-of-magnitude luminosities 
of collapsing objects, assuming a central temperature T ~ 1000 K. 

In table |] we list the luminosities of 5 important lines, at all the four stages we tabulated 
in table [Bl~| , plus the total line luminosity, the continuum luminosity, and their sum. They 
are compared to the values predicted by KS01. 

It can be seen that KS01 predict a total luminosity comparable to our results for the 
initial phases of collapse; instead, in the last phases our total luminosities are £ 10 times 
larger than the KS01 value. This is quite reasonable, since the approximations they used 
are much better at the beginning of the collapse rather than at the end. If we consider the 
luminosities of single lines, a more important discrepancy emerges. KS01 correctly predict 
that the rotational transition at A = 9.69/im is the most important one at the beginning, but 
they estimate a very high line luminosity, that no lines actually achieve, even in advanced 
stages. The reason is that they attribute too large a fraction of the total luminosity to single 
lines. Instead, the emission is actually split into a growing number of lines: in the early 
phases, the two rotational lines listed in table [| accounts for ~ 50% of total emission (first 
tabulated stage), but the luminosities of single lines "saturate" when their central optical 
depth becomes £ 1. As the total line emission rises, we have mainly an increase in the 
number of important lines, rather than a growth of the luminosity of single lines: at the last 
stage no line accounts for more than ~ 5% of total line emission. 

In table ^ we also give the expected specific fluxes coming from a single protostellar 
object at a redshift z = 19 in an SCDM cosmological model at the first and the, for the first 
and the last listed stage. It can be seen that the most promising lines are those produced 
by purely rotational transitions, such as those at A = 9.69 /im and A = 8.06 /im (which are 
redshifted to A G bs — 190 and A Q bs — 160 /im, respectively), whose expected specific flux is 
£ 10 nJy for at least ~ 10 4 yr); vibrational transitions, such as those at A = 2.41 /im and 
A = 2.81 /xhi (redshifted to A b s — 48 /im and A b s — 56 /im) reach specific fluxes of the same 
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Table 4. Comparison of line and continuum luminosities; luminosities (L) are in units of 10 34 ergs - 1 , specific fluxes at earth 
(1, for an object at z = 19 in a SCDM cosmology with Hq = 75 kms — 1 Mpc — x ) in njy; the continuum specific flux is calculated 
at the peak of the black-body spectrum. (A ~ 2fim). All data come from run P100/2 
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- 


0.00 


0.13 


75 


- 


- 


0.002 


Total 


4.4 


43.2 


89.3 


181 


6.7 


- 


- 



order only in the last phases and their emission at such levels should last only for £, 10 3 yr. 
Finally, we predict that the continuum specific flux should be of the order of 0.002 njy at an 
observed wavelength of X b s ~ 40yum (from a rest-frame wavelength of ~ 2/im); we cannot 
predict the duration of continuum emission, since it is deeply affected by the evolution of 
the accretion and by the evolution of the hydrostatic core at the centre. 

Such fluxes are largely below the detection threshold of any observational facility that 
will be available in the near future. NGST sensitivity for continuum radiation should be of 
the order of ~ 1 /iJy at A = 20 /zm, while ALMA should reach a sensitivity of ~ 10 mJy for 
a line 25 kms -1 wide at A ~ 350 /xm. However, the possibility that primordial protostars 
form in large groups must not to be excluded, since a galactic-scale (10 11 M & ) fluctuation 
could host a large number of 10 6 M Q halos, and each of them could host several protostellar 
objects. For instance KS01 roughly estimate that a galactic-scale fluctuation could have a 
luminosity ~ 10 4 times larger than a single object, and that such luminosity could last for 
~ 10 7 yr. If this is correct, specific fluxes from galactic-scale objects should be of the order 
of ~ 0.01 yuJy and ~ 0.1 mJy for continuum and line emission, respectively. 

5 CONCLUSIONS 

In this work we have investigated the evolution of the luminosity and spectra of the radiation 
emitted by primordial collapsing protostars. We find that initially radiation comes mostly 
from H2 rotational lines. Later, the contributions from H2 vibrational lines and also from 
continuum radiation (due to the accretion upon a centrally hydrostatic core) become impor- 
tant. At early stages, only two rotational lines account for about 50% of the total emission, 
but the subsequent increase in total line luminosity results in an increase of the number 
of prominent lines, rather than a large increase in the luminosity of the most important 
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ones. As a result, the expected line specific fluxes are always ^ 50 nJy (with line widths 
£ 30 kms -1 ), even for the most important lines, which should fall at observed wavelengths 
in the range 150 — 200/im. Continuum radiation should reach a specific flux of the order of 
0.002 nJ, peaking at an observed wavelength between 20 and 40 /im. These fluxes should be 
quite universal, provided that a molecular central region of ~ 3M develops at the centre 
of each protostellar cloud. The predicted fluxes, both for line and continuum radiation, are 
largely below the detection threshold of existing or planned instruments (such as NGST and 
ALMA), but the possibility that primordial protostars form large "clusters" remains to be 
investigated, even if rough estimates do not look very promising (KS01). 

In our work we also found that, different from the present-day case, radiation forces play 
virtually no role in the initial evolution of a protostellar object (i.e. , during the self-similar 
collapse, the formation of the core and the onset of accretion). Although we were not able 
to completely bridge the gap with accretion models such as those of SPS86 and Omukai & 
Palla 2001, we have shown that their initial assumptions are quite reasonable. Even if there 
are differences in the ionization structure, there is indication (in particular, the decrease in 
the mass accretion rate) that these will vanish as the hydrostatic core mass grows. Thus, 
undisturbed spherical collapse should proceed until the initial protostellar cloud has been 
largely accreted, or until radiation forces halt the infall at M* ~ 300 M Q (Omukai & Palla 
2001). 

What can stop mass growth is the presence of HD molecules (even in a tiny fraction) 
that may result from shocks (Flower & Pineau des Forets 2001, Uehara & Inutsuka 2000); 
in this case the fragmentation of the PopIII object might end up in clumps of roughly solar 
rather than thousand solar masses, thus drastically changing the final mass of the star. 

If this process is not important and angular momentum and magnetic fields do not play 
a significant role either, the final mass of a primordial protostar is primarily set by the 
amount of gas available in the parent cloud. Then, it is conceivable that the primordial IMF 
is skewed towards high masses, and that very massive stars (VMS), with masses in excess of 
M ~ 100 Mq, can form. Such stars may end-up as a population of intermediate mass black 
holes (Portinari et al. 1998, Heger & Langer 2000) that represent the end product of such 
an episode of pregalactic star formation. Then, these seeds could be incorporated into the 
larger units as seen in the bulges of present-day galaxies (Madau & Rees 2001). 

We thank an anonymous referee for a critical reading and of the manuscript and for 
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stimulating remarks. We also thank B. Ciardi for discussion and for making her chemistry 
routines available to us, and K. Omukai, R. Nishi and F. Palla for discussion and clarifying 
remarks about their works. 

APPENDIX A: THE CODE 

Hereafter we describe the main features of the code we developed. 

Al Generalities 

A 1.1 Gas hydrodynamics 

Despite the cosmological context of our calculations, the collapsing baryon clouds we are 
studying are assumed to be self-gravitating, and we can neglect the gravitational effects of 
host dark matter halo. 

The code integrates the following 1-D Lagrangian hydrodynamical equations (conserva- 
tion of mass, momentum and energy) using an explicit finite-difference scheme: 

dM r A 2 , . 1 . 

— — = Anry (Al) 

dr 

dv n dp GM r , . . . 

Tt = - 4 "' lk - — + f " d (A2) 

du d ( l\ V — A ... 

* = '"a [p) + —< (A3) 

where the Lagrangian coordinate M r is the mass inside a radius r. p, v, p, f ra d, u, represent 
the density, velocity, pressure, radiative force per unit mass and internal energy per unit 
mass, respectively. T and A are the heating and cooling rate per unit volume, respectively. 
Hydrodynamical effects are included by adding to the pressure an artificial viscosity term 
of the form 

qocpAv 2 , (A4) 

in regions where contiguous gas shells would compenetrate (e.g. Bowers & Wilson 1991; 
Thoul & Weinberg 1995; hereafter, TW95). 

A 1.2 Chemistry 

Because of the extremely low metallicity of the primordial material as predicted by standard 
big-bang nucleosynthesis (Z ~ 10~ 12 ; Kolb & Turner 1990), we can safely limit our treat- 
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Table Al. Chemical reactions and reaction rates 
Reaction Rate 



H+ + e" -» 


H + 7 


AANZ98, k 2 


H + e- -> 


H" +7 


AANZ98, k 7 


H + H" -» 


H 2 + e- 


AANZ98, fc g 


H + H + H -> 


H 2 +H 


PSS83, fc 4 


H 2 +H^ 


H + H + H 


PSS83, fc 5 


H + H + H 2 — 


H 2 +H 2 


PSS83, fc 6 


H 2 + H 2 -> 


H + H + H 2 


PSS83, fc 7 


H + e- -> 


H+ + c - + e _ 


AANZ98, fci 


H + H ->• 


H+ + H + e- 


PSS83, k 9 



AANZ98 stands for Abel, Anninos, Norman & Zhang (1998) and PSS83 stands for Palla, Salpeter & Stahlcr (1983). Note that 
the H recombination rate (first reaction in the list) is appropriate for case A (case B could be more appropriate, but results 
are not affected by the difference, as can be seen from the comparison to ON98). 

ment of chemistry to the various ionic species of atomic hydrogen (H°, H + , H~), molecular 
hydrogen (H2 , Hj), and helium (He , He + , He ++ ), plus free electrons (e~). 

The He number abundance is taken to be 1/12 of the total H number abundance; this 
translates into mass fractions of X=0.75, Y=0.25 and Z=0. 

Some authors (Flower 2000; Flower & Pineau de Forets 2001; and also, in a different 
context, Uehara & Inutsuka 2000) suggested that deuterium is an important coolant in a 
metal-free gas, because of the permanent dipole moment of HD molecules and their low 
excitation temperature (~ 160 K for HD instead of ~ 510 K for H 2 ). However HD cooling 
is likely to be important only at temperatures below some hundred degrees, i.e. during 
the formation of protostellar clouds, rather than during their collapse. For this reason, we 
neglect the effects of deuterium. Note that BCL01 included HD cooling in their investigation 
of primordial cloud fragmentation, and found that (for a standard Deuterium abundance) 
the thermal evolution of the gas would not change significatively 

Initial conditions were chosen with chemical abundances far from equilibrium, and we 
follow their evolution by integrating a network of reactions (see Table 1). In particular, we 
include the 3-body reactions described by Palla, Salpeter & Stahler (1983). 

As the gas gets denser (p £ p cq = 5 x 10~ n g cm -3 ), the abundances of the relevant 
species (H2 , H° and He , which then represent more than 99.9% of the total number abun- 
dance) are close to their equilibrium values, as reactions. We then switch to the equilibrium 
chemistry. 

When the gas is almost ideal, i.e. when p < 10 4 dyne cm -2 (p is the gas pressure) 
we find the equilibrium abundances by solving the Saha equations for the relevant species. 
However these equations become inappropriate at high temperatures and densities, so when 
p > 10 6 dyne cm -2 we interpolate the abundance values from the tables given by SCVH95. 
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In the intermediate regime (10 4 dyne cm -3 < p < 10 6 dyne cm -2 ), we evaluate the chemical 
composition both with Saha equations and using the SCVH95 tables, and then we appropri- 
ately interpolate the two results; this is done in order to ensure a smooth transition between 
the two regimes. 

A 1.3 Thermodynamical properties 

For the solution of the hydrodynamical equations, we need to specify the equation of state 
and to calculate the internal energy as a function of temperature, density and chemical 
composition. 

At low densities and temperatures (in practice, we use p < 10 4 dyne cm -2 as a criterion), 
the gas is ideal to a very good level of approximation and we can use the usual equation 

p = nkT (A5) 

where n is the number density of the gas. 

In this ideal regime we can also write the total internal energy per unit mass as 

U = Mth + Mchem (A6) 

with 

«th = —^—N tot kT, (A7) 

7ad - 1 

9 

Mchcm = -N tot Y^ ViXi (A8) 

where u t h and Wchcm are the thermal and the chemical energy per unit mass, while A tot , yi 
and Xi are the total number of particles per unit mass, the number fraction of particles of 
species i (y, = Ni/N tot ), and the binding energy of a particle of species i, respectively. Last, 
7 a d is the mean adiabatic index of the gas, obtained by a weighted average of the adiabatic 
indexes of the various chemical species: for the seven included mono-atomic species we have 
7ad,x = 5/3, while for the two diatomic species we use 7 ad H + = 7/5 and the expression (from 
ON98) is 

1 5 2 e x 6100K /A „ X 

= - + x 2 - -7T, x = . (A9) 

lad,H 2 -l 2 (e--l) 2 ' T l ; 

The value of 7 for H^ neglects the vibrational degrees of freedom, but the abundance of this 
ion is always very small and this does not introduce any significant error. 

When the gas is no more ideal (p > 10 6 dyne cm -2 ), we find p(p,T) and u(p,T) by 
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interpolating the values tabulated by SCVH95, which have been calculated jeeping non- 
ideal effects into account. In the transition between the two regimes, we interpolate the 
results of the ideal and of the SCVH95 values in order to obtain a smooth transition. 

A 1.4 Radiative transfer and cooling 

In the initial stages of collapse, gas cooling is dominated by the H 2 roto-vibrational transi- 
tions, which soon become optically thick. A careful treatment of radiative transfer effects is 
then crucial in the study of the collapse. 

H 2 line cooling is important only when T £ 2500 K, as at higher temperatures H 2 molecules 
slowly dissociate, and continuum emission starts dominating the local cooling. The energy 
gaps between H2 rotational levels typically correspond to excitation temperatures of fewx 100 
K - 10 3 K, while the gap between vibrational levels corresponds to ~ 6000 K. Under such 
conditions, it is adequate to consider only H 2 roto-vibrational levels with rotational quantum 
number < J < 20, and vibrational quantum number < v < 2. In this way, all the states 
whose excitation temperatures (relative to the ground state) is ^ 20000 K are considered. 
We treat the radiative transfer of all the lines (between the considered levels) permitted by 
the quadrupole selection rules. 

Line opacities are computed using the Einstein A-coefhcients from Turner, Kirby-Docken 
& Dalgarno (1977), and assuming local thermodynamic equilibrium (LTE) for estimating 
the level populations. This is correct since the density at which a line gives a relevant 
contribution to the total cooling is always greater than the critical density at which collisional 
de-excitations dominate upon radiative transitions. 

Finally, the effects of thermal broadening and Doppler shifts due to bulk motions arising 
in the collapse phase have been included in the line transfer calculations. 

When the temperature rises to a sufficiently high value (T > 1000 K), we include the 
cooling by continuum radiation; continuum radiative transfer is treated as a gray problem, 
using the mean Planck opacity for Z = gas from Lenzuni, Chernoff & Salpeter (1991) when 
1000 K < T < 7000 K, and the Rosseland mean opacity for metal-free gas from Rogers & 
Iglesias (1992) when T > 7000 K. 

A2 Numerical details 
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A2.1 Implementation of the basic equations 

The numerical treatment of the dynamical evolution is made along the lines described by 
Bowers & Wilson (1991), and TW95. 

With the exception of the energy equation, the dynamical equations listed in §2.1 were 
discretized as in TW95 (their equations 9, 10, 11, 12, 13 and 15), assuming a value of 4 for 
the artificial viscosity parameter c q . 

Such Lagrangian finite difference scheme is second order accurate in space and time. 

For the energy equation (|A3|) , our treatment starts from equation (14) of TW95. The 
chemical term in the internal energy (see §2.3) introduces a complex dependence upon the 
final temperature. For this reason the energy equation is solved iteratively, using the bisection 
method. In this case the result is only first order accurate in time. 

A2.2 Time-steps 

The integration time-step is chosen evaluating the minimum among: 

- twice the previous integration time-step; 

- the minimum dynamical time-scale over all shells i (mmi{Cd yn (Grf / M r>i y/ 2 }); 

- the minimum Courant time-scale over all shells i (mmi{Ccom^ r i/ Cs,i})', 

- the minimum cooling time-scale over all shells i (mmi{C coo \Uip/ (Ti — Aj)}); 

- the minimum crossing time-scale over all shells i (mmi{C CIOSS Ari / Avi}); 

In addition, we check a posteriori that the temperature does not change by more than 0.5% 
in a single time-step. The factors Cd yn = 1, Cc OU r = 0.2, C coo i = 0.1 and C cross = 0.05 are 
safety constants. 

A2.3 Boundary conditions 

As internal boundary condition we simply require that the centre of our object is a geomet- 
rical point (r = 0, M r = M = 0), with zero velocity (t> in = 0). 

The external boundary condition is set by taking a constant pressure outside the most 
external shell. 

A2.4 Implementation of radiative transfer 
General scheme 
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The solution of the radiative transfer equation in spherical symmetry requires to find the 
specific intensity I v (r,/j,), where /i is the cosine of the angle 9 between the radial outward 
direction and the direction of the considered ray. 

This is done by means of the standard tangent ray method (Hummer & Rybicki 1971; 
Bowers & Wilson 1991; ON98). In this scheme, the specific intensity I u (r, fi) is evaluated at 
the external radius rj of each shell for 1 + 2i values of \i. 

Once I(ri,fij) is known, we compute the monochromatic flux and the net outward bolo- 
metric luminosity L(ri) at each shell radius. 

The cooling rate (per unit time per unit volume) of the shell i is then given by 
_ L{n) - Irfo-Q 
(4vr/3)(r.f -rti) 

H 2 lines 

To calculate the source function and the absorption coefficient, we assume LTE. This is 
a fair assumption, as the density is always much greater than the relevant critical values. 
The absorption coefficient is given by 

hu , s {2hv 3 „ \ / g 2 \ . k , 

<*» = 4^") (-T^i J (nx| - n 2 j (All) 

where 4>{y) is the line profile, A 2 \ is the Einstein coefficient for the relevant transition (taken 
from Turner et al. 1977), g 2 and g\ are the statistical weights of the two levels involved in 
the transition, and n 2 and n x are number densities of molecules in the two considered levels, 
computed assuming LTE. 

As noted by ON98, since the temperature typical of H2 line emission is quite low, the line 
profiles are narrow, and the Doppler shift caused by the gas bulk infall can be important. 
We assume a Gaussian profile for each line, sampled in 30 frequency bins, and we properly 
take into account Doppler shift in the line transfer. If the line absorption depth is < 0.1, we 
use the optically thin approximation. 
Continuum 

The continuum transfer is treated as a gray problem. Again, we assume LTE to estimate 
the source function. The absorption coefficient is calculated from the (Planck) opacity tables 
of Lenzuni et al. 1991 when T < 7000 K, and from the (Rosseland) opacity tables of Rogers 
& Iglesias 1992 for T > 7000 K. 

As the continuum optical depth can reach values 3> 1 (see fig. J8|), the numerical treatment 
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above can become inaccurate. For this reason, when r > 4 in a single shell (see the conditions 
listed in ON98), we use the diffusion approximation, i.e. 
F - ^aT*dT^_ maTfT^-TU 
3a dr 3«j r i+1 - n-i 

where a is the Stefan-Boltzmann constant. 

APPENDIX B: LIST OF EMITTED LINES 

In the following table we list the luminosity at infinity of each H2 line we included into our 
code. Lines are sorted according to their wavelength. 
The columns give, respectively: 

- the quantum levels describing the transition; 

- the wavelength (in /iin) of the line; 

- the radial optical depth from the centre of the simulated object to infinity, as evaluated 
at line centre, and at the epoch when the continumm optical depth is r c = 1; 

- the line total luminosity at infinity, evaluated at four different epochs, i.e. when the 
central temperature is 650 K (L65o)> 1000 K (L 1000 ) and 1500 K (L 1500 ), and when r c = 1 
(L Tc= i); such luminosities are given in erg s -1 ; 

- the specific luminosity (in erg s" 1 Hz -1 ) at the stage when r c = 1, estimated using a 
line width equal to the Doppler broadening Au = (u /c)(2kT/ /imp}) 1 / 2 , where u is the line 
central frequency, fi ~ 2 is the H 2 molecular weight and T is the temperature of the shell 
where the radial optical depth to infinity at line centre is equal to 1 (or T ~ 2338 K for the 
few lines with r < 1). 

The last lines of the table give the total line luminosities, the continuum luminosity and 
the sum of the two. Continuum luminosity is not calculated (being actually very small) 
when the central temperature is < 1000 K; continuum specific luminosity is only a rough 
estimate, since we took Az/ = f pea k, where z/ pca k — 1.4 x 10 14 Hz is the frequency of the peak 
of blackbody emission at T = 2338 K. 
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Table Bl. Line list 



(VI, Ji) 


-> (Vf, J{) 


A[/im] 


T 


^650 


^1000 


^1500 


^T C 


=1 


'r c =l 


(2,11) - 


-(0,9) 


1.0700 


4.39 x 10 1 


9.91 


x 


10 24 


1.59 x 


10 29 


6.50 x 10 31 


4.55 x 


10 32 


9.77 x 10 22 


(2 


10)- 


-(0,8) 


1.0701 


7.65 x 10 1 


7.00 


x 


10 25 


6.32 x 


10 29 


1.54 x 10 32 


6.50 x 


10 32 


1.40 x 10 23 


(2 


12)- 


- (0, 10) 


1.0736 


2.34 x 10 1 


1.24 


x 


10 24 


3.64 x 


10 28 


2.36 x 10 31 


2.85 x 


10 32 


6.13 x 10 22 


(2 


9)- 


'(0,7) 


1.0737 


1.24 x 10 2 


4.32 


X 


10 26 


2.25 x 


10 30 


2.97 x 10 32 


9.08 x 


10 32 


2.26 x 10 23 


(2 


8)- 


■(0,6) 


1.0806 


1.83 x 10 2 


2.25 


X 


10 27 


6.99 x 


10 30 


4.84 x 10 32 


1.19 x 


10 33 


2.97 x 10 23 


(2 


13)- 


-(o,n) 


1.0809 


1.18 x 10 1 


1.43 


X 


10 23 


7.73 x 


10 27 


8.01 x 10 30 


1.76 x 


10 32 


3.48 x 10 22 


(2 


7)- 


.(0,5) 


1.0908 


2.46 x 10 2 


9.68 


X 


10 27 


1.85 x 


10 31 


6.81 x 10 32 


1.44 x 


10 33 


3.65 x 10 23 


(2 


14)- 


- (0, 12) 


1.0925 


5.57 x 10° 


1.55 


X 


10 22 


1.54 x 


10 27 


2.58 x lO 30 


1.05 x 


10 32 


2.09 x 10 22 


(2 


6)- 


>(0,4) 


1.1042 


2.98 x 10 2 


3.36 


X 


10 28 


3.91 x 


10 31 


8.61 x 10 32 


1.67 x 


10 33 


4.26 x 10 23 


(2 


15)- 


- (0, 13) 


1.1087 


2.51 x 10° 


1.61 


X 


10 21 


2.97 x 


10 26 


7.98 x 10 29 


5.79 x 


10 31 


1.18 x 10 22 


(2 


5)~ 


.(0,3) 


1.1210 


3.19 x 10 2 


9.17 


X 


10 28 


6.66 x 


10 31 


1.02 x 10 33 


1.89 x 


10 33 


4.91 x 10 23 


(2 


16)- 


- (0, 14) 


1.1300 


1.07 x 10° 


1.65 


X 


10 20 


5.53 x 


10 25 


2.38 x 10 29 


2.80 x 


10 31 


5.79 x 10 21 


(2 


4)- 


'(0,2) 


1.1410 


2.96 x 10 2 


1.90 


X 


10 29 


9.34 x 


10 31 


1.19 x 10 33 


2.06 x 


10 33 


5.45 x 10 23 


(2 


17)- 


- (0, 15) 


1.1573 


0.43 x 10° 


1.69 


X 


10 19 


1.00 x 


10 25 


6.79 x 10 28 


1.19 x 


10 31 


2.52 x 10 21 


(2 


3)- 


'(0,1) 


1.1645 


2.27 x 10 2 


2.81 


X 


10 29 


1.09 x 


10 32 


1.29 x 10 33 


2.20 x 


10 33 


5.93 x 10 23 


(2 


2)- 


.(0,0) 


1.1914 


1.25 x 10 2 


2.59 


X 


10 29 


9.28 x 


10 31 


1.24 x 10 33 


2.13 x 


10 33 


5.88 x 10 23 


(2 


18)- 


- (0, 16) 


1.1915 


0.16 x 10° 


1.78 


X 


10 18 


1.81 x 


10 24 


1.89 x 10 28 


4.54 x 


10 30 


9.91 x 10 2() 


(2 


19)- 


- (0, 17) 


1.2338 


5.57 x 10~ 2 


1.88 


X 


10 17 


3.14 x 


10 23 


4.87 x 10 27 


1.55 x 


10 30 


3.49 x 10 20 


(2 


1)- 


'(0,1) 


1.2400 


1.20 x 10 2 


3.59 


X 


10 29 


1.17 x 


10 32 


1.47 x 10 33 


2.44 X 


10 33 


7.01 x 10 23 


(2 


2)- 


'(0,2) 


1.2439 


1.23 x 10 2 


2.69 


X 


10 29 


1.02 x 


10 32 


1.44 x 10 33 


2.46 X 


10 33 


7.08 x 10 23 


(2 


3)- 


'(0,3) 


1.2498 


1.31 x 10 2 


1.78 


X 


10 29 


8.42 x 


10 31 


1.37 x 10 33 


2.40 X 


10 33 


6.94 x 10 23 


(2 


4)- 


'(0,4) 


1.2577 


1.24 x 10 2 


9.04 


X 


10 28 


5.69 x 


10 31 


1.23 x 10 33 


2.24 X 


10 33 


6.53 x 10 23 


(2 


5)- 


'(0,5) 


1.2676 


1.06 x 10 2 


3.60 


X 


10 28 


3.13 x 


10 31 


1.00 x 10 33 


2.01 X 


10 33 


5.89 x 10 23 


(2 


6)- 


.(0,6) 


1.2796 


8.25 x 10 1 


1.14 


X 


10 28 


1.45 x 


10 31 


7.18 x 10 32 


1.67 x 


10 33 


4.30 x 10 23 


(2 


20) - 


- (0, 18) 


1.2860 


1.58 x 10~ 2 


1.88 


X 


10 16 


4.89 x 


10 22 


1.08 x 10 27 


4.39 x 


10 29 


1.03 x 10 20 


(2 


7)- 


'(0,7) 


1.2938 


6.00 x 10 1 


2.99 


X 


10 27 


5.78 x 


10 30 


4.51 x 10 32 


1.30 x 


10 33 


3.37 x 10 23 


(2 


0)- 


'(0,2) 


1.2947 


6.98 x 10 1 


2.58 


X 


10 29 


8.42 x 


10 31 


1.30 x 10 33 


2.28 x 


10 33 


6.82 x 10 23 


(2 


8)- 


.(0,8) 


1.3101 


4.04 x 10 1 


6.45 


X 


10 26 


2.01 x 


10 30 


2.38 x 10 32 


9.45 x 


10 32 


2.49 x 10 23 


(2 


9)- 


.(0,9) 


1.3287 


2.59 x 10 1 


1.20 


X 


10 26 


6.24 x 


10 29 


1.10 x 10 32 


6.46 x 


10 32 


1.72 x 10 23 


(2 


1)- 


'(0,3) 


1.3372 


8.54 x 10 1 


2.76 


X 


10 29 


9.68 x 


10 31 


1.52 x 10 33 


2.65 x 


10 33 


8.19 x 10 23 


(2 


10)- 


- (0, 10) 


1.3497 


1.57 x 10 1 


1.94 


X 


10 25 


1.75 x 


10 29 


4.64 x 10 31 


4.30 x 


10 32 


1.16 x 10 23 


(2 


11)- 


-(0,11) 


1.3730 


9.12 x 10° 


2.80 


X 


10 24 


4.50 x 


10 28 


1.85 x 10 31 


2.68 x 


10 32 


6.74 x 10 22 


(2 


2)~ 


.(0,4) 


1.3838 


7.50 x 10 1 


1.81 


X 


10 29 


7.26 x 


10 31 


1.42 x 10 33 


2.52 x 


10 33 


8.08 x 10 23 


(2 


12)- 


- (0, 12) 


1.3988 


5.14 x 10° 


3.72 


X 


10 23 


1.09 x 


10 28 


7.05 x 10 30 


1.70 x 


10 32 


4.35 x 10 22 


(2 


13)- 


- (0, 13) 


1.4271 


2.81 x 10° 


4.64 


X 


10 22 


2.50 x 


10 27 


2.59 x 10 30 


9.89 x 


10 31 


2.58 x 10 22 


(2 


3)- 


'(0,5) 


1.4349 


5.51 x 10 1 


8.39 


X 


10 28 


4.11 x 


10 31 


1.08 x 10 33 


2.13 x 


10 33 


6.14 x 10 23 


(2 


14)- 


- (0, 14) 


1.4581 


1.51 x 10° 


5.63 


X 


10 21 


5.62 x 


10 26 


9.38 x 10 29 


5.56 x 


10 31 


1.48 x 10 22 


(2 


4)- 


'(0,6) 


1.4903 


3.49 x 10 1 


2.88 


X 


10 28 


1.82 x 


10 31 


6.56 x 10 32 


1.55 x 


10 33 


4.65 x 10 23 


(2 


15)- 


- (0, 15) 


1.4918 


0.81 x 10° 


6.88 


X 


10 20 


1.26 x 


10 26 


3.40 x 10 29 


3.01 x 


10 31 


8.22 x 10 21 


(2 


16)- 


- (0, 16) 


1.5283 


0.44 x 10° 


8.69 


X 


10 19 


2.90 x 


10 25 


1.25 x 10 29 


1.59 x 


10 31 


4.46 x 10 21 


(2 


5)- 


'(0,7) 


1.5503 


1.97 x 10 1 


7.54 


X 


10 27 


6.56 x 


10 30 


3.12 x 10 32 


1.01 x 


10 33 


3.14 x 10 23 


(2 


17)- 


- (0, 17) 


1.5676 


0.24 x 10° 


1.17 


X 


10 19 


6.96 x 


10 24 


4.72 x 10 28 


8.39 x 


10 30 


2.41 x 10 21 


(2 


18)- 


- (0, 18) 


1.6097 


0.14 x 10° 


1.76 


X 


10 18 


1.79 x 


10 24 


1.86 x 10 28 


4.49 x 


10 30 


1.32 x 10 21 


(2 


6)- 


'(0,8) 


1.6146 


9.92 x 10° 


1.54 


X 


10 27 


1.94 x 


10 30 


1.21 x 10 32 


5.81 x 


10 32 


1.88 x 10 23 


(2 


19)- 


- (0, 19) 


1.6546 


7.92 x 10~ 2 


3.03 


X 


10 17 


5.08 x 


10 23 


7.86 x 10 27 


2.50 x 


10 30 


7.57 x 10 2() 


(2 


7)- 


'(0,9) 


1.6830 


4.58 x 10° 


2.51 


X 


10 26 


4.85 x 


10 29 


4.12 x 10 31 


2.90 x 


10 32 


8.93 x 10 22 


(1 


13)- 


-(o,n) 


1.6989 


9.37 x 10° 


2.45 


X 


10 25 


1.66 x 


10 29 


3.55 x 10 31 


3.58 x 


10 32 


1.11 x 10 23 


(1 


12)- 


- (0, 10) 


1.7019 


3.36 x 10 1 


4.43 


X 


10 26 


1.54 x 


10 30 


2.01 x 10 32 


9.96 x 


10 32 


3.40 x 10 23 


(2 


20) - 


- (0, 20) 


1.7022 


4.86 x 10~ 2 


6.29 


X 


10 16 


1.64 x 


10 23 


3.60 x 10 27 


1.47 x 


10 30 


4.57 x 10 2() 


(1 


14)- 


- (0, 12) 


1.7048 


1.58 x 10° 


8.06 


X 


10 23 


1.07 x 


10 28 


3.79 x 10 30 


8.03 x 


10 31 


2.51 x 10 22 


(1 


11)- 


-(0,9) 


1.7134 


9.49 x 10 1 


6.07 


X 


10 27 


1.10 x 


10 31 


8.47 x 10 32 


2.21 x 


10 33 


7.61 x 10 23 


(1 


15)- 


- (0, 13) 


1.7203 


2.19 x 10~ 2 


2.19 


X 


10 21 


5.67 x 


10 25 


3.35 x 10 28 


1.12 x 


10 30 


3.52 x 10 20 


(1 


10)- 


-(0,8) 


1.7331 


2.25 x 10 2 


6.57 


X 


10 28 


6.38 x 


10 31 


2.43 x 10 33 


4.32 x 


10 33 


1.73 x 10 24 


(1 


16)- 


- (0, 14) 


1.7464 


0.23 x 10° 


4.68 


X 


10 21 


2.34 x 


10 26 


2.28 x 10 29 


1.07 x 


10 31 


3.42 x 10 21 


(2 


8)- 


' (0, 10) 


1.7552 


1.92 x 10° 


3.29 


X 


10 25 


1.02 x 


10 29 


1.23 x 10 31 


1.27 x 


10 32 


4.08 x 10 22 


(1 


9)- 


'(0,7) 


1.7610 


4.66 x 10 2 


5.76 


X 


10 29 


3.09 x 


10 32 


5.01 x 10 33 


7.29 x 


10 33 


2.97 x 10 24 


(1 


17)- 


- (0, 15) 


1.7847 


0.65 x 10° 


2.76 


X 


10 21 


2.60 x 


10 26 


4.15 x 10 29 


2.62 x 


10 31 


8.55 x 10 21 


(1 


8)- 


'(0,6) 


1.7974 


8.55 x 10 2 


4.05 


X 


10 30 


1.16 x 


10 33 


8.32 x 10 33 


1.13 x 


10 34 


4.72 x 10 24 


(2 


12)- 


-(1,10) 


1.8248 


2.55 x 10° 


1.85 


X 


10 23 


5.41 x 


10 27 


3.51 x 10 30 


1.02 x 


10 32 


3.39 x 10 22 


(2 


13)- 


-(1,11) 


1.8259 


0.40 x 10° 


6.58 


X 


10 21 


3.55 x 


10 26 


3.68 x 10 29 


1.77 x 


10 31 


5.93 x 10 21 


(2 


9)- 


.(0,11) 


1.8302 


0.75 x 10° 


3.60 


X 


10 24 


1.87 x 


10 28 


3.30 x 10 30 


4.74 x 


10 31 


1.59 x 10 22 


(2 


11)- 


-(1,9) 


1.8333 


8.88 x 10° 


2.76 


X 


10 24 


4.44 x 


10 28 


1.82 x 10 31 


2.81 x 


10 32 


9.42 x 10 22 


(2 


14)- 


-(1,12) 


1.8371 


1.74 x 10~ 4 


6.35 


X 


10 17 


6.33 x 


10 22 


1.06 x 10 26 


7.16 x 


10 27 


2.41 x 10 18 


(1 


18)- 


- (0, 16) 


1.8372 


0.89 x 10° 


8.63 


X 


10 20 


1.48 x 


10 26 


3.77 x 10 29 


3.20 x 


10 31 


1.08 x 10 22 
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Table Bl continued 



(Wi.Ji)- 


-> (Vf, J{) 


A[/im] 


T 




^650 


^1000 


^1500 


Lt c 


=1 


'r c =l 


(1,7)- 


(0,5) 


1.8426 


1.39 x 


10 3 


2.26 


X 


IO 31 


3.04 x 


IO 33 


1.26 x IO 34 


1.55 x 


IO 34 


8.12 x IO 24 


(2 


10)- 


-(1,8) 


1.8511 


2.34 x 


10 1 


2.96 


x 


IO 25 


2.67 x 


IO 29 


7.08 x IO 31 


6.21 x 


IO 32 


2.31 x IO 23 


(2 


15)- 


-(1,13) 


1.8595 


0.13 x 


10° 


1.03 


X 


IO 20 


1.90 x 


IO 25 


5.11 x IO 28 


4.80 x 


IO 30 


1.64 x IO 21 


(2 


9)- 


(1,7) 


1.8780 


5.15 x 


10 1 


2.45 


X 


IO 26 


1.28 x 


IO 30 


2.25 x IO 32 


1.18 x 


IO 33 


4.44 x IO 23 


(2 


16)- 


- (1, 14) 


1.8944 


0.31 x 


10° 


5.73 


X 


IO 19 


1.91 x 


IO 25 


8.24 x IO 28 


1.06 x 


IO 31 


3.69 x IO 21 


(1 


6)- 


(0,4) 


1.8971 


1.99 x 


10 3 


9.77 


X 


IO 31 


5.83 x 


IO 33 


1.68 x IO 34 


2.01 x 


IO 34 


1.09 x IO 25 


(1 


19)- 


- (0, 17) 


1.9070 


0.97 x 


10° 


2.37 


X 


IO 20 


7.08 x 


IO 25 


2.79 x IO 29 


3.16 x 


IO 31 


1.10 x IO 22 


(2 


10)- 


- (0, 12) 


1.9071 


0.27 x 


10° 


3.32 


X 


IO 23 


3.00 x 


IO 27 


7.96 x IO 29 


1.55 x 


IO 31 


5.42 x IO 21 


(2 


8)- 


(1,6) 


1.9142 


9.77 x 


10 1 


1.61 


X 


IO 27 


5.00 x 


IO 30 


5.94 x IO 32 


2.02 x 


IO 33 


7.78 x IO 23 


(2 


17)- 


-(1,15) 


1.9439 


0.41 x 


10° 


1.85 


X 


IO 19 


1.10 x 


IO 25 


7.45 x IO 28 


1.31 x 


IO 31 


4.65 x IO 21 


(2 


7)- 


(1,5) 


1.9600 


1.62 x 


10 2 


8.29 


X 


IO 27 


1.61 x 


IO 31 


1.27 x IO 33 


3.06 x 


IO 33 


1.20 x IO 24 


(1 


5)- 


(0,3) 


1.9617 


2.51 x 


10 3 


3.21 


X 


IO 32 


9.35 x 


IO 33 


2.20 x IO 34 


2.55 x 


IO 34 


1.43 x IO 25 


(2 


11)- 


- (0, 13) 


1.9845 


8.94 x 


io- 2 


2.65 


X 


IO 22 


4.27 x 


IO 26 


1.75 x IO 29 


4.59 x 


IO 30 


1.67 x IO 21 


(1 


20) - 


- (0, 18) 


1.9985 


0.95 x 


10° 


6.70 


X 


IO 19 


3.30 x 


IO 25 


1.93 x IO 29 


2.87 x 


IO 31 


1.05 x IO 22 


(2 


18)- 


-(1,16) 


2.0108 


0.44 x 


10° 


5.19 


X 


10 1S 


5.29 x 


IO 24 


5.51 x IO 28 


1.30 x 


IO 31 


4.79 x IO 21 


(2 


6)- 


(1,4) 


2.0160 


2.35 x 


10 2 


3.34 


X 


IO 28 


4.23 x 


IO 31 


2.21 x IO 33 


4.37 x 


IO 33 


2.04 x IO 24 


(1 


4)- 


(0,2) 


2.0372 


2.70 x 


10 3 


7.62 


X 


IO 32 


1.28 x 


IO 34 


2.68 x IO 34 


2.95 x 


IO 34 


1.71 x IO 25 


(2 


12)- 


- (0, 14) 


2.0602 


2.80 x 


10~ 2 


1.88 


X 


IO 21 


5.50 x 


IO 25 


3.57 x IO 28 


1.27 x 


IO 30 


4.78 x IO 20 


(2 


5)- 


(1,3) 


2.0829 


2.98 x 


10 2 


1.03 


X 


IO 29 


8.95 x 


IO 31 


3.17 x IO 33 


5.33 x 


IO 33 


2.57 x IO 24 


(2 


19)- 


-(1,17) 


2.0994 


0.43 x 


10° 


1.45 


X 


IO 18 


2.43 x 


IO 24 


3.76 x IO 28 


1.17 x 


IO 31 


4.51 x IO 21 


(1 


3)- 


(0,1) 


2.1249 


2.40 x 


10 3 


1.28 


X 


IO 33 


1.62 x 


IO 34 


3.05 x IO 34 


3.40 x 


IO 34 


2.06 x IO 25 


(2 


13)- 


- (0, 15) 


2.1320 


8.21 x 


L0~ 3 


1.21 


X 


IO 20 


6.52 x 


IO 24 


6.76 x IO 27 


3.31 x 


IO 29 


1.29 x IO 20 


(2 


4)-, 


(1,2) 


2.1617 


3.25 x 


10 2 


2.36 


X 


IO 29 


1.49 x 


IO 32 


3.97 x IO 33 


6.30 x 


IO 33 


3.15 x IO 24 


(2 


14)- 


- (0, 16) 


2.1968 


2.27 x 


10~ 3 


7.26 


X 


IO 18 


7.25 x 


IO 23 


1.21 x IO 27 


8.20 x 


IO 28 


3.30 x IO 19 


(2 


20) - 


-(1,18) 


2.2155 


0.41 x 


10° 


4.45 


X 


IO 17 


1.16 x 


IO 24 


2.55 x IO 28 


1.02 x 


IO 31 


4.13 x IO 21 


(1 


2)- 


(0,0) 


2.2261 


1.54 x 


10 3 


1.33 


X 


IO 33 


1.67 x 


IO 34 


3.19 x IO 34 


3.53 x 


IO 34 


2.24 x IO 25 


(2 


15)- 


- (0, 17) 


2.2514 


6.02 x 


io- 4 


4.22 


X 


IO 17 


7.76 x 


IO 22 


2.09 x IO 26 


1.96 x 


IO 28 


8.10 x IO 18 


(2 


3)- 


(1,1) 


2.2537 


2.91 x 


10 2 


3.81 


X 


IO 29 


1.86 x 


IO 32 


4.20 x IO 33 


6.44 x 


IO 33 


3.36 x IO 24 


(2 


16)- 


- (0, 18) 


2.2922 


1.55 x 


io- 4 


2.48 


X 


IO 16 


8.29 x 


IO 21 


3.57 x IO 25 


4.64 x 


IO 27 


1.95 x IO 18 


(2 


17)- 


- (0, 19) 


2.3159 


4.12 x 


10~ 5 


1.60 


X 


IO 15 


9.48 x 


IO 20 


6.42 x IO 24 


1.15 x 


IO 27 


4.86 x IO 17 


(2 


18)- 


- (0, 20) 


2.3196 


1.26 x 


IO" 5 


1.29 


X 


IO 14 


1.32 x 


IO 20 


1.37 x IO 24 


3.31 x 


IO 26 


1.40 x IO 17 


(2 


2)-. 


(1,0) 


2.3603 


1.88 x 


IO 2 


3.78 


X 


IO 29 


1.52 x 


IO 32 


3.38 x IO 33 


5.33 x 


IO 33 


2.91 x IO 24 


(1 


1)- 


(0,1) 


2.4094 


1.82 x 


IO 3 


2.04 


X 


IO 33 


2.21 x 


IO 34 


3.90 x IO 34 


4.26 x 


IO 34 


2.92 x IO 25 


(1 


2)- 


(0,2) 


2.4165 


1.85 x 


IO 3 


1.49 


X 


IO 33 


2.04 x 


IO 34 


3.64 x IO 34 


3.98 x 


IO 34 


2.74 x IO 25 


(1 


3)- 


(0,3) 


2.4273 


1.92 x 


IO 3 


9.28 


X 


IO 32 


1.77 x 


IO 34 


3.41 x IO 34 


3.76 x 


IO 34 


2.60 x IO 25 


(1 


4)-. 


(0,4) 


2.4416 


1.78 x 


IO 3 


4.34 


X 


IO 32 


1.36 x 


IO 34 


3.03 x IO 34 


3.27 x 


IO 34 


2.27 x IO 25 


(1 


5)- 


(0,5) 


2.4597 


1.49 x 


IO 3 


1.56 


X 


IO 32 


9.01 x 


IO 33 


2.37 x IO 34 


2.72 x 


IO 34 


1.91 x IO 25 


(1 


6)- 


(0,6) 


2.4815 


1.13 x 


IO 3 


4.37 


X 


IO 31 


4.65 x 


IO 33 


1.81 x IO 34 


2.08 x 


IO 34 


1.20 x IO 25 


(1 


7)- 


(0,7) 


2.5071 


7.88 x 


IO 2 


9.81 


X 


IO 30 


1.81 x 


IO 33 


1.19 x IO 34 


1.46 x 


IO 34 


8.49 x IO 24 


(1 


8)- 


(0,8) 


2.5366 


5.12 x 


IO 2 


1.81 


X 


IO 30 


5.59 x 


IO 32 


7.21 x IO 33 


9.57 x 


IO 33 


5.62 x IO 24 


(2 


1)- 


(1,1) 


2.5547 


2.26 x 


IO 2 


5.72 


X 


IO 29 


2.02 x 


IO 32 


3.92 x IO 33 


5.77 x 


IO 33 


3.41 x IO 24 


(2 


2)- 


(1,2) 


2.5633 


2.34 x 


IO 2 


4.26 


X 


IO 29 


1.71 x 


IO 32 


3.74 x IO 33 


5.77 x 


IO 33 


3.42 x IO 24 


(1 


9)- 


(0,9) 


2.5701 


3.13 x 


IO 2 


2.81 


X 


IO 29 


1.51 x 


IO 32 


3.84 x IO 33 


5.84 x 


IO 33 


3.47 x IO 24 


(2 


3)- 


(1,3) 


2.5763 


2.46 x 


IO 2 


2.76 


X 


IO 29 


1.35 x 


IO 32 


3.41 x IO 33 


5.40 x 


IO 33 


3.22 x IO 24 


(2 


4)-, 


(1,4) 


2.5937 


2.32 x 


IO 2 


1.37 


X 


IO 29 


8.68 x 


IO 31 


2.80 x IO 33 


4.59 x 


IO 33 


2.76 x IO 24 


(1 


10)- 


- (0, 10) 


2.6077 


1.80 x 


IO 2 


3.72 


X 


IO 28 


3.60 x 


IO 31 


1.71 x IO 33 


3.49 x 


IO 33 


1.83 x IO 24 


(2 


5)- 


(1,5) 


2.6156 


1.98 x 


IO 2 


5.31 


X 


IO 28 


4.62 x 


IO 31 


1.99 x IO 33 


3.72 x 


IO 33 


1.95 x IO 24 


(1 


0)- 


(0,2) 


2.6297 


1.32 x 


IO 3 


1.63 


X 


IO 33 


2.18 x 


IO 34 


3.79 x IO 34 


4.11 x 


IO 34 


3.08 x IO 25 


(2 


6)- 


(1,6) 


2.6421 


1.54 x 


IO 2 


1.64 


X 


IO 28 


2.07 x 


IO 31 


1.24 x IO 33 


2.90 x 


IO 33 


1.54 x IO 24 


(1 


11)- 


-(o,n) 


2.6495 


9.92 x 


IO 1 


4.37 


X 


IO 27 


7.92 x 


IO 30 


6.37 x IO 32 


1.91 x 


IO 33 


1.02 x IO 24 


(2 


7)- 


(1,7) 


2.6735 


1.11 X 


IO 2 


4.09 


X 


IO 27 


7.93 x 


IO 30 


6.68 x IO 32 


1.99 x 


IO 33 


1.07 x IO 24 


(1 


12)- 


- (0, 12) 


2.6955 


5.26 x 


IO 1 


4.64 


X 


IO 26 


1.61 x 


IO 30 


2.11 x IO 32 


9.93 x 


IO 32 


5.37 x IO 23 


(2 


8)- 


(1,8) 


2.7097 


7.41 x 


IO 1 


8.51 


X 


IO 26 


2.65 x 


IO 30 


3.18 x IO 32 


1.28 x 


IO 33 


6.94 x IO 23 


(1 


13)- 


- (0, 13) 


2.7459 


2.71 x 


IO 1 


4.60 


X 


IO 25 


3.11 x 


IO 29 


6.67 x IO 31 


5.20 x 


IO 32 


2.87 x IO 23 


(2 


9)- 


(1,9) 


2.7510 


4.68 x 


IO 1 


1.50 


X 


IO 26 


7.81 x 


IO 29 


1.38 x IO 32 


8.02 x 


IO 32 


4.43 x IO 23 


(2 


0)- 


(1,2) 


2.7891 


1.69 x 


IO 2 


4.45 


X 


IO 29 


1.47 x 


IO 32 


2.99 x IO 33 


4.57 x 


IO 33 


2.95 x IO 24 


(2 


10)- 


-(1,10) 


2.7976 


2.81 x 


IO 1 


2.31 


X 


IO 25 


2.08 x 


IO 29 


5.52 x IO 31 


4.79 x 


IO 32 


2.69 x IO 23 


(1 


14)- 


- (0, 14) 


2.8008 


1.36 x 


IO 1 


4.37 


X 


IO 24 


5.79 x 


IO 28 


2.06 x IO 31 


2.71 x 


IO 32 


1.39 x IO 23 


(1 


1)- 


(0,3) 


2.8057 


1.87 x 


IO 3 


1.77 


X 


IO 33 


2.49 x 


IO 34 


4.22 x IO 34 


4.55 x 


IO 34 


3.64 x IO 25 


(2 


11)- 


-(1,11) 


2.8498 


1.60 x 


IO 1 


3.14 


X 


IO 24 


5.04 x 


IO 28 


2.07 x IO 31 


2.76 x 


IO 32 


1.44 x IO 23 


(1 


15)- 


- (0, 15) 


2.8601 


6.79 x 


10° 


4.14 


X 


IO 23 


1.07 x 


IO 28 


6.34 x IO 30 


1.40 x 


IO 32 


7.35 x IO 22 


(2 


12)- 


-(1,12) 


2.9077 


8.82 x 


10° 


3.90 


X 


IO 23 


1.14 x 


IO 28 


7.39 x IO 30 


1.62 x 


IO 32 


8.63 x IO 22 


(1 


16)- 


- (0, 16) 


2.9240 


3.39 x 


10° 


4.03 


X 


IO 22 


2.01 x 


IO 27 


1.97 x IO 30 


7.09 x 


IO 31 


3.80 x IO 22 


(2 


13)- 


- (1, 13) 


2.9716 


4.73 x 


10° 


4.56 


X 


IO 22 


2.46 x 


IO 27 


2.55 x IO 30 


8.82 x 


IO 31 


4.80 x IO 22 
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(v 


,Ji)- 


-> (v f ,Jf) 


A [/»m] 


T 


^650 


^1000 


^1500 


Lr c 


=1 


'r c =l 


(2 


1)- 


(1,3) 


2.9774 


2.44 X 10 2 


4.94 


X 


IO 29 


1.74 x 


IO 32 


3.51 x 


IO 33 


5.26 x 


IO 33 


3.62 x IO 24 


(1 


17)- 


- (0, 17) 


2.9923 


1.70 x 10" 


4.16 


x 


IO 21 


3.92 x 


IO 26 


6.26 x 


IO 29 


3.57 x 


IO 31 


1.96 x IO 22 


(1 


2)- 


(0,4) 


3.0073 


1.92 x IO 3 


1.16 


x 


IO 33 


2.24 X 


IO 34 


3.94 x 


IO 34 


4.24 x 


IO 34 


3.64 x IO 25 


(2 


14)- 


-(1,14) 


3.0417 


2.50 x 10" 


5.17 


x 


IO 21 


5.16 x 


IO 26 


8.61 x 


IO 29 


4.72 x 


IO 31 


2.63 x IO 22 


(1 


18)- 


- (0, 18) 


3.0651 


0.88 x 10" 


4.78 


x 


IO 20 


8.20 x 


IO 25 


2.09 x 


IO 29 


1.78 x 


IO 31 


10.00 x IO 21 


(2 


15)- 


-(1,15) 


3.1184 


1.30 x 10" 


5.80 


x 


IO 20 


1.07 x 


IO 26 


2.87 x 


IO 29 


2.43 x 


IO 31 


1.39 x IO 22 


(1 


19)- 


•+ (0, 19) 


3.1422 


0.47 x 10" 


6.28 


x 


IO 19 


1.88 x 


IO 25 


7.40 x 


IO 28 


8.86 x 


IO 30 


5.10 x IO 21 


(2 


2)- 


(1,4) 


3.1941 


2.56 x 10 2 


3.35 


x 


IO 29 


1.35 x 


IO 32 


3.07 x 


IO 33 


4.63 x 


IO 33 


3.42 x IO 24 


(2 


16)- 


-(1,16) 


3.2018 


0.68 x 10" 


6.72 


x 


IO 19 


2.25 x 


IO 25 


9.66 x 


IO 28 


1.20 x 


IO 31 


7.05 x IO 21 


(1 


20) - 


- (0, 20) 


3.2232 


0.26 x 10" 


9.88 


x 


IO 18 


4.86 x 


IO 24 


2.84 x 


IO 28 


4.55 x 


IO 30 


2.69 x IO 21 


(1 


3)- 


(0,5) 


3.2385 


1.65 x 10 3 


5.24 


x 


IO 32 


1.61 x 


IO 34 


3.22 x 


IO 34 


3.37 x 


IO 34 


3.11 x IO 25 


(2 


17)- 


-(1,17) 


3.2923 


0.36 x 10" 


8.30 


x 


IO 18 


4.92 x 


IO 24 


3.33 x 


IO 28 


5.88 x 


IO 30 


3.54 x IO 21 


(2 


18)- 


-(1,18) 


3.3899 


0.19 x 10" 


1.12 


x 


IO 18 


1.14 x 


IO 24 


1.19 x 


IO 28 


2.86 x 


IO 30 


1.78 x IO 21 


(2 


3)-» 


(1,5) 


3.4438 


2.25 x 10 2 


1.59 


x 


IO 29 


7.79 x 


IO 31 


2.18 x 


IO 33 


3.57 x 


IO 33 


2.47 x IO 24 


(2 


19)- 


-(1,19) 


3.4948 


0.11 x 10" 


1.75 


x 


IO 17 


2.93 x 


IO 23 


4.53 x 


IO 27 


1.44 x 


IO 30 


9.22 x IO 2 " 


(1 


4)-. 


(0,6) 


3.5034 


1.22 x 10 3 


1.71 


x 


IO 32 


8.63 x 


IO 33 


2.22 x 


IO 34 


2.37 x 


IO 34 


1.92 x IO 25 


(2 


20) - 


-(1,20) 


3.6070 


6.16 x 10~ 2 


3.23 


x 


IO 16 


8.39 x 


IO 22 


1.85 x 


IO 27 


7.53 x 


IO 29 


4.97 x IO 2 " 


(2 


4)- 


(1,6) 


3.7313 


1.72 x 10 2 


5.56 


x 


IO 28 


3.52 x 


IO 31 


1.26 x 


IO 33 


2.37 x 


IO 33 


1.78 x IO 24 


(1 


5)- 


(0,7) 


3.8063 


8.07 x IO 2 


4.18 


x 


IO 31 


3.29 x 


IO 33 


1.33 x 


IO 34 


1.49 x 


IO 34 


1.31 x IO 25 


(2 


5)- 


(1,7) 


4.0622 


1.18 x IO 2 


1.48 


x 


IO 28 


1.29 x 


IO 31 


6.06 x 


IO 32 


1.45 x 


IO 33 


1.19 x IO 24 


(1 


6)- 


(0,8) 


4.1517 


4.84 x IO 2 


7.87 


x 


IO 30 


9.33 x 


IO 32 


6.68 x 


IO 33 


7.78 x 


IO 33 


7.48 x IO 24 


(° 


15)- 


- (0, 13) 


4.3167 


1.04 x IO 3 


5.78 


X 


IO 27 


1.71 x 


IO 31 


1.63 x 


IO 33 


3.08 x 


IO 33 


3.08 x IO 24 


(° 


16)- 


- (0, 14) 


4.3364 


6.68 x IO 2 


6.49 


x 


IO 26 


3.93 x 


IO 30 


7.34 x 


IO 32 


1.86 x 


IO 33 


1.62 x IO 24 


(° 


14)- 


- (0, 12) 


4.3568 


1.59 x IO 3 


5.08 


X 


IO 28 


7.30 x 


IO 31 


3.22 x 


IO 33 


4.75 x 


IO 33 


4.79 x IO 24 


(° 


17)- 


- (0, 15) 


4.4219 


4.35 x IO 2 


7.52 


X 


IO 25 


9.09 x 


IO 29 


3.00 x 


IO 32 


1.18 x 


IO 33 


1.05 x IO 24 


(2 


6)- 


(1,8) 


4.4422 


7.29 x IO 1 


3.05 


X 


IO 27 


3.86 x 


IO 30 


2.40 x 


IO 32 


7.47 x 


IO 32 


6.66 x IO 23 


(° 


13)- 


- (0,11) 


4.4558 


2.38 x IO 3 


4.26 


X 


IO 29 


2.96 x 


IO 32 


5.79 x 


IO 33 


7.21 x 


IO 33 


7.44 x IO 24 


(1 


7)~> 


(0,9) 


4.5431 


2.63 x IO 2 


1.15 


X 


IO 30 


2.15 x 


IO 32 


2.80 x 


IO 33 


3.74 x 


IO 33 


3.93 x IO 24 


(° 


18)- 


- (0, 16) 


4.5859 


2.89 x IO 2 


9.34 


X 


IO 24 


2.18 x 


IO 29 


1.19 x 


IO 32 


7.07 x 


IO 32 


6.51 x IO 23 


(1 


15)- 


- (1,13) 


4.6058 


1.07 x IO 2 


2.65 


X 


IO 24 


6.86 x 


IO 28 


4.05 x 


IO 31 


3.80 x 


IO 32 


3.52 x IO 23 


(° 


12)- 


- (0, 10) 


4.6173 


3.43 x IO 3 


3.27 


X 


IO 3 " 


1.11 X 


IO 33 


9.48 x 


IO 33 


1.08 x 


IO 34 


1.16 x IO 25 


(1 


14)- 


- (1,12) 


4.6386 


1.58 x IO 2 


1.97 


X 


IO 25 


2.62 x 


IO 29 


9.29 x 


IO 31 


5.62 x 


IO 32 


5.23 x IO 23 


(1 


16)- 


- (1,14) 


4.6391 


7.20 x IO 1 


3.55 


X 


IO 23 


1.77 x 


IO 28 


1.73 x 


IO 31 


2.43 x 


IO 32 


2.27 x IO 23 


(1 


13)- 


-(1,11) 


4.7356 


2.28 x IO 2 


1.40 


X 


IO 26 


9.44 x 


IO 29 


2.02 x 


IO 32 


8.39 x 


IO 32 


7.98 x IO 23 


(1 


17)- 


- (1, 15) 


4.7462 


4.87 x IO 1 


4.86 


X 


IO 22 


4.58 x 


IO 27 


7.31 x 


IO 30 


1.61 x 


IO 32 


1.53 x IO 23 


(° 


11)- 


- (0, 9) 


4.8495 


4.68 x IO 3 


2.20 


X 


IO 31 


3.64 x 


IO 33 


1.49 x 


IO 34 


1.65 x 


IO 34 


1.86 x IO 25 


(° 


19) - 


- (0, 17) 


4.8513 


2.02 x IO 2 


1.30 


X 


IO 24 


5.58 x 


IO 28 


4.86 x 


IO 31 


4.45 x 


IO 32 


4.33 x IO 23 


(2 


7)- 


(1,9) 


4.8764 


4.10 x IO 1 


4.95 


X 


IO 26 


9.59 x 


IO 29 


8.15 x 


IO 31 


3.54 x 


IO 32 


3.46 x IO 23 


(1 


12)- 


- (i, io) 


4.8999 


3.17 x IO 2 


9.11 


X 


IO 26 


3.17 x 


IO 30 


4.07 x 


IO 32 


1.13 x 


IO 33 


1.11 x IO 24 


(1 


18)- 


- (1, 16) 


4.9429 


3.40 x IO 1 


7.16 


X 


IO 21 


1.23 x 


IO 27 


3.13 x 


IO 30 


1.03 x 


IO 32 


9.32 x IO 22 


(2 


15)- 


- (2, 13) 


4.9684 


1.31 x IO 1 


2.31 


X 


IO 21 


4.24 x 


IO 26 


1.14 x 


IO 30 


5.49 x 


IO 31 


4.99 x IO 22 


(1 


8)- 


(0,10) 


4.9824 


1.31 x IO 2 


1.35 


X 


IO 29 


4.17 x 


IO 31 


9.06 x 


IO 32 


1.52 x 


IO 33 


1.52 x IO 24 


(2 


14)- 


- (2, 12) 


4.9894 


1.86 x IO 1 


1.44 


X 


IO 22 


1.44 x 


IO 27 


2.40 x 


IO 30 


7.54 x 


IO 31 


6.89 x IO 22 


(2 


16)- 


- (2, 14) 


5.0221 


9.20 x 10" 


3.67 


X 


IO 20 


1.23 x 


IO 26 


5.27 x 


IO 29 


3.84 x 


IO 31 


3.53 x IO 22 


(2 


13)- 


- (2,H) 


5.0817 


2.57 x IO 1 


8.60 


X 


IO 22 


4.64 x 


IO 27 


4.80 x 


IO 30 


1.02 x 


IO 32 


9.49 x IO 22 


(1 


11)- 


-(1,9) 


5.1400 


4.19 x IO 2 


5.27 


X 


IO 27 


9.57 x 


IO 30 


7.26 x 


IO 32 


1.52 x 


IO 33 


1.57 x IO 24 


(2 


17)- 


- (2, 15) 


5.1610 


6.57 x 10" 


6.04 


X 


IO 19 


3.59 x 


IO 25 


2.43 x 


IO 29 


2.64 x 


IO 31 


2.49 x IO 22 


(° 


10)- 


- (0, 8) 


5.1673 


5.98 x IO 3 


1.26 


X 


IO 32 


9.17 x 


IO 33 


2.09 x 


IO 34 


2.29 x 


IO 34 


3.38 x IO 25 


(2 


12)- 


- (2, 10) 


5.2478 


3.45 x IO 1 


4.75 


X 


IO 23 


1.39 x 


IO 28 


9.01 x 


IO 30 


1.32 x 


IO 32 


1.40 x IO 23 


(1 


19)- 


- (1, 17) 


5.2578 


2.45 x IO 1 


1.16 


X 


IO 21 


3.46 x 


IO 26 


1.36 x 


IO 30 


6.60 x 


IO 31 


6.36 x IO 22 


(° 


20) - 


- (0, 18) 


5.2594 


1.51 x IO 2 


2.09 


X 


IO 23 


1.56 x 


IO 28 


2.09 x 


IO 31 


2.70 x 


IO 32 


2.85 x IO 23 


(2 


8)- 


(1,10) 


5.3686 


2.12 x IO 1 


6.51 


X 


IO 25 


2.02 x 


IO 29 


2.44 x 


IO 31 


1.56 x 


IO 32 


1.69 x IO 23 


(2 


18)- 


- (2, 16) 


5.4058 


4.82 x 10" 


1.06 


X 


IO 19 


1.08 x 


IO 25 


1.13 x 


IO 29 


1.79 x 


IO 31 


1.77 x IO 22 


(1 


9)- 


(0,11) 


5.4679 


6.01 x IO 1 


1.29 


X 


IO 28 


6.91 x 


IO 30 


2.32 x 


IO 32 


5.77 x 


IO 32 


6.34 x IO 23 


(1 


10)- 


-(1,8) 


5.4711 


5.17 x IO 2 


2.60 


X 


IO 28 


2.52 x 


IO 31 


1.12 x 


IO 33 


1.90 x 


IO 33 


2.40 x IO 24 


(2 


11)- 


- (2, 9) 


5.4959 


4.40 x IO 1 


2.34 


X 


IO 24 


3.77 x 


IO 28 


1.55 x 


IO 31 


1.56 x 


IO 32 


1.73 x IO 23 


(° 


9)- 


(0,7) 


5.5942 


6.98 x IO 3 


5.82 


X 


IO 32 


1.72 x 


IO 34 


2.80 x 


IO 34 


2.97 x 


IO 34 


4.73 x IO 25 


(1 


20) - 


- (1, 18) 


5.7428 


1.96 x IO 1 


2.22 


X 


IO 20 


1.09 x 


IO 26 


6.38 x 


IO 29 


4.37 x 


IO 31 


4.60 x IO 22 


(2 


19)- 


- (2, 17) 


5.7939 


3.71 x 10" 


2.06 


X 


IO 18 


3.45 x 


IO 24 


5.34 x 


IO 28 


1.19 x 


IO 31 


1.27 x IO 22 


(2 


10)- 


- (2, 8) 


5.8420 


5.26 x IO 1 


9.91 


X 


IO 24 


8.94 x 


IO 28 


2.37 x 


IO 31 


1.86 x 


IO 32 


2.18 x IO 23 


(1 


9)- 


(1,7) 


5.9180 


5.88 x IO 2 


1.05 


X 


IO 29 


5.61 x 


IO 31 


1.51 x 


IO 33 


2.17 x 


IO 33 


2.97 x IO 24 


(2 


9)-> 


(1,11) 


5.9189 


1.01 x IO 1 


7.05 


X 


IO 24 


3.67 x 


IO 28 


6.46 x 


IO 30 


6.33 x 


IO 31 


6.86 x IO 22 


(1 


10)- 


- (0, 12) 


5.9915 


2.55 x IO 1 


1.03 


X 


IO 27 


1.00 x 


IO 30 


5.15 x 


IO 31 


1.98 x 


IO 32 


2.38 x IO 23 


(° 


8)- 


(0,6) 


6.1682 


7.36 x IO 3 


2.09 


X 


IO 33 


2.61 x 


IO 34 


3.58 x 


IO 34 


3.52 x 


IO 34 


6.19 x IO 25 
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Table Bl continued 



(vi,Ji) -* (vt, J f ) A [/mi] 



-^650 



ir c = l 



l Tc = l 



(2.9) - 


►(2,7) 


6.3119 


5.83 x 10 1 


3.48 


x 10 25 


1.81 


x 


10 29 


3.19 


X 


10 31 


1.92 


X 


10 32 


2.44 


x 10 23 


(2, 20) - 


- (2, 18) 


6.3949 


3.09 x 10° 


4.59 


x 10 17 


1.19 


X 


10 24 


2.63 


X 


10 28 


7.87 


X 


10 30 


9.21 


xlO 21 


(1,8)- 


►(1.6) 


6.5205 


6.06 x 10 2 


3.31 


x 10 29 


1.02 


X 


10 32 


1.72 


X 


10 33 


2.23 


X 


10 33 


3.36 


xlO 24 


(2, 10) - 


-(1,12) 


6.5207 


4.46 x 10° 


6.44 


x 10 23 


5.81 


X 


10 27 


1.54 


X 


10 30 


2.33 


X 


10 31 


2.78 


x 10 22 


(1,11)- 


-* (0, 13) 


6.5354 


9.87 x 10° 


7.05 


x 10 25 


1.28 


X 


10 29 


1.04 


X 


10 31 


6.46 


X 


10 31 


7.73 


x 10 22 


(2.8) - 


,(2,6) 


6.9478 


5.83 x 10 1 


9.62 


x 10 25 


2.99 


X 


10 29 


3.60 


X 


10 31 


1.74 


X 


10 32 


2.43 


x 10 23 


(0, 7) - 


.(0,5) 


6.9523 


6.86 x 10 3 


5.52 


x 10 33 


3.30 


X 


10 34 


3.87 


X 


10 34 


3.92 


X 


10 34 


7.77 


x 10 25 


(1,12)- 


- (0, 14) 


7.0689 


3.51 x 10° 


4.23 


xlO 24 


1.47 


X 


10 28 


1.92 


X 


10 30 


1.95 


X 


10 31 


2.53 


x 10 22 


(2,11)- 


-(1,13) 


7.1563 


1.81 x 10° 


5.08 


x 10 22 


8.17 


X 


10 26 


3.36 


X 


10 29 


7.78 


X 


10 30 


1.02 


x 10 22 


(1,7)- 


,(1,5) 


7.3448 


5.49 x 10 2 


7.79 


xlO 29 


1.45 


X 


10 32 


1.66 


X 


10 33 


2.06 


X 


10 33 


3.51 


xlO 24 


(1,18)- 


- (0, 20) 


7.3467 


4.62 x 10" 4 


3.78 


x 10 16 


6.48 


X 


10 21 


1.65 


X 


10 25 


1.50 


X 


10 27 


2.01 


x 10 18 


(1,13) - 


- (0, 15) 


7.5460 


1.13 x 10° 


2.28 


x 10 23 


1.54 


X 


10 27 


3.30 


X 


10 29 


5.28 


X 


10 30 


7.30 


x 10 21 


(2, 12) - 


-(1,14) 


7.7921 


0.67 x 10° 


3.51 


x 10 21 


1.03 


X 


10 26 


6.66 


X 


10 28 


2.29 


X 


10 30 


3.27 


x 10 21 


(1,17) - 


- (0, 19) 


7.8090 


3.27 x 10" 3 


9.91 


x 10 17 


9.33 


X 


10 22 


1.49 


X 


10 26 


9.79 


X 


10 27 


1.40 


xlO 19 


(2, 7) - 


,(2,5) 


7.8200 


5.17 x 10 1 


2.01 


x 10 26 


3.89 


X 


10 29 


3.31 


X 


10 31 


1.39 


X 


10 32 


2.18 


x 10 23 


(1,14) - 


- (0, 16) 


7.9088 


0.33 x 10° 


1.14 


x 10 22 


1.51 


X 


10 26 


5.35 


X 


10 28 


1.26 


X 


10 30 


1.83 


x 10 21 


(0, 6) - 


,(0,4) 


8.0563 


5.44 x 10 3 


9.89 


x 10 33 


3.48 


X 


10 34 


3.72 


X 


10 34 


3.74 


X 


10 34 


8.58 


x 10 25 


(1,16)- 


- (0, 18) 


8.0684 


1.81 x 10" 2 


2.37 


xlO 19 


1.18 


X 


10 24 


1.16 


X 


10 27 


5.43 


X 


10 28 


8.02 


xlO 19 


(1,15)- 


- (0, 17) 


8.0977 


8.27 x 10" 2 


5.32 


xlO 20 


1.38 


X 


10 25 


8.13 


X 


10 27 


2.72 


X 


10 29 


4.03 


xlO 20 


(2, 18) - 


-(1,20) 


8.2737 


8.43 x 10~ 5 


6.52 


x 10 13 


6.64 


X 


10 19 


6.91 


X 


10 23 


1.67 


X 


10 26 


2.52 


x 10 17 


(2,13) - 


-(1,15) 


8.3749 


0.22 x 10° 


2.20 


x 10 20 


1.18 


X 


10 25 


1.23 


X 


10 28 


5.98 


X 


10 29 


9.18 


x 10 20 


(1,6)- 


►(1,4) 


8.5070 


4.27 x 10 2 


1.28 


x 10 30 


1.52 


X 


10 32 


1.31 


X 


10 33 


1.56 


X 


10 33 


3.08 


xlO 24 


(2, 17) - 


-(1,19) 


8.8069 


6.51 x 10~ 4 


1.61 


x 10 15 


9.55 


X 


10 20 


6.47 


X 


10 24 


1.15 


X 


10 27 


1.86 


x 10 18 


(2, 14) - 


-(1,16) 


8.8339 


6.61 x 10~ 2 


1.27 


x 10 19 


1.26 


X 


10 24 


2.11 


X 


10 27 


1.43 


X 


10 29 


2.31 


x 10 20 


(2,6) - 


►(2,4) 


9.0514 


3.94 x 10 1 


2.97 


x 10 26 


3.75 


X 


10 29 


2.34 


X 


10 31 


9.27 


X 


10 31 


1.68 


x 10 23 


(2, 16) - 


-(1,18) 


9.0901 


3.72 x 10" 3 


3.45 


x 10 16 


1.15 


X 


10 22 


4.96 


X 


10 25 


6.45 


X 


10 27 


1.07 


x 10 19 


(2,15) - 


-(1,17) 


9.0924 


1.71 x 10~ 2 


6.82 


x 10 17 


1.25 


X 


10 23 


3.37 


X 


10 26 


3.17 


X 


10 28 


5.29 


x 10 19 


(0, 5) - 


,(0,3) 


9.6894 


3.54 x 10 3 


1.12 


x 10 34 


2.78 


X 


10 34 


2.90 


X 


10 34 


2.90 


X 


10 34 


7.99 


x 10 25 


(1,5)- 


►(1,3) 


10.2273 


2.73 x 10 2 


1.35 


x 10 30 


1.08 


X 


10 32 


7.61 


X 


10 32 


8.85 


X 


10 32 


2.10 


x 10 24 


(2,5)- 


,(2,3) 


10.8762 


2.47 x 10 1 


2.86 


x 10 26 


2.49 


X 


10 29 


1.18 


X 


10 31 


4.56 


X 


10 31 


9.96 


x 10 22 


(0,4)- 


,(0,2) 


12.3006 


1.75 x 10 3 


7.30 


x 10 33 


1.58 


X 


10 34 


1.55 


X 


10 34 


1.53 


X 


10 34 


5.35 


x 10 25 


(1,4)- 


►(1,2) 


12.9795 


1.32 x 10 2 


7.97 


x 10 29 


4.57 


X 


10 31 


2.93 


X 


10 32 


3.39 


X 


10 32 


1.02 


xlO 24 


(2,4) - 


,(2,2) 


13.7974 


1.19 x 10 1 


1.56 


x 10 26 


9.89 


X 


10 28 


3.75 


X 


10 30 


1.46 


X 


10 31 


3.69 


x 10 22 


(0, 3) - 


► (0, 1) 


17.0595 


5.63 x 10 2 


2.08 


x 10 33 


4.73 


X 


10 33 


4.58 


X 


10 33 


4.47 


X 


10 33 


2.17 


x 10 25 


(1,3)- 


►(1,1) 


17.9969 


4.22 x 10 1 


2.03 


xlO 29 


8.87 


X 


10 30 


5.23 


X 


10 31 


6.75 


X 


10 31 


2.44 


x 10 23 


(2,3) - 


► (2, 1) 


19.1253 


3.73 x 10° 


3.72 


x 10 25 


1.82 


X 


10 28 


5.73 


X 


10 29 


2.49 


X 


10 30 


8.72 


x 10 21 


(0, 2 ) - 


► (0, 0) 


28.2559 


7.90 x 10 1 


1.27 


x 10 32 


3.51 


X 


10 32 


3.63 


X 


10 32 


3.36 


X 


10 32 


2.20 


xlO 24 


(1,2)- 


►(1,0) 


29.8042 


5.88 x 10° 


1.12 


x 10 28 


4.01 


X 


10 29 


2.10 


X 


10 30 


3.35 


X 


10 30 


1.83 


xlO 22 


(2,2) - 


,(2,0) 


31.6666 


0.52 x 10° 


1.96 


xlO 24 


7.89 


X 


10 26 


2.16 


X 


10 28 


9.82 


X 


10 28 


5.69 


x 10 20 


sum of lines 


— 


— 


5.32 


xlO 34 


4.32 


X 


10 35 


8.92 


X 


10 35 


1.06 


X 


10 36 




— 


continuum 


— 


1.00 x 10° 




— 


8.39 


X 


10 27 


1.29 


X 


10 33 


7.50 


X 


10 35 


~ 5 


x 10 21 


total 




— 


— 


5.32 


xlO 34 


4.32 


X 


10 35 


8.93 


X 


10 35 


1.81 


X 


10 36 




— 
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